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Executive  Summary 


The  research  objective  for  Phase  I  of  this  work  was  to  demonstrate  the  feasibility  of 
establishing  optimal  shapes  for  hypervelocity  projectiles,  considering  simultaneously  the 
effects  of  drag,  internal  volume  distribution  and  impact  effectiveness.  Utilizing  classical 
optimization  analysis  and  hypersonic  inviscid  flow  theory,  minimum  pressure  drag  bodies 
were  obtained  subject  to  various  constraints.  It  was  then  assumed  that  the  shapes  so 
obtained  corresponded  to  the  solid  body  plus  a  viscous  displacement.  Utilizing  state-of- 
the-art  turbulent  boundary  layer  theory  in  an  iterative  mode,  the  solid  body  shape  was 
then  determined  along  with  the  skin  friction  drag.  Relative  to  any  specified  reference 
penetrator  shape,  the  internal  volume  distribution  effectiveness  may  then  be  character¬ 
ized.  Subsequently,  the  loss  in  kinetic  energy  of  the  projectile  along  its  trajectory  Is 
computed  and  the  initial  stage  of  Impact  is  computed  from  a  Navier-Stokes  numerical 
simulation  In  order  to  quantify  the  target  damage.  As  a  result  of  this  procedure,  various 
candidate  hypervelocity  projectile  shapes  may  be  competed  against  one  another  from  an 
overall  mission  point  of  view- 


Discussion 


The  prediction  of  optimum  aerodynamic  shapes  has  a  heritage  dating  from  Newton  to 
the  present  time.  Traditionally,  simplified  inviscid  flow  models  are  combined  with  classi¬ 
cal  optimization  theory  in  order  to  derive  minimum  drag  configurations1.  More  recently, 
incorporation  of  viscous  effects  at  on-  and  off-design  conditions  has  been  approached, 
as  shown  in  the  Appendix.  Relative  to  the  smart  hypervelocity  projectile  design  problem, 
additional  considerations  beyond  aerodynamic  drag  are  relevant.  Included  in  these  are 
the  internal  volume  distribution  and  the  impact  mechanics  of  the  penetrator  which  fits 
within  that  volume.  The  present  research  sought  to  make  a  contribution  to  the  design 
of  optimally-shaped  hypervelocity  projectiles  by  considering  all  aspects  of  the  mission, 
including  viscous  effects  on  the  minimum  drag  body,  as  shown  in  Fig.  1. 

Initially,  the  minimum  pressure  dreg  shape  (herein  assumed  to  correspond  to  the  solid 
body  plus  its  viscous  displacement)  for  various  constraints  may  be  sought  from  standard 
optimization  techniques  applied  to  the  integral: 

1  =  /  cp  d» 

where  Cp  Is  the  local  pressure  coefficient  and  c&  is  the  local  surface  area  is  the  stream- 
wise  (drag)  direction  divided  by  the  base  area.  Kcco.ding  to  the  Eul*r-Lagrange  equation, 
the  pressure  drag  is  minimized  when: 


d/d*  (30/3r‘)  -  3j/3r  =  o 


» 


J  =  Cp  dA/dx 
dA/dx  =  r  r' 

T  =  r/rbase 
r ' =  dr/dx 

For  axisymmetrlc  flow,  the  Newton-Busemarm  approximation  for  the  local  pressure  in 
hypersonic  flow  may  be  written: 

Cp  *  2.0  (0*  +  r0*/2) 

where  ©  Is  the  local  inclination  of  the  surface  to  the  flow.  This  simple  formula  allows 
the  Euler  >  Lagrange  equation  to  be  integrated  subject  to  various  constraints,  as  we 
shall  see.  In  addition,  however,  it  is  possible  to  integrate  the  Euler  -  Lagrange  equation 
numerically  and  compute  the  local  pressure  field  from  hypersonic  small  disturbance  theory, 
wherein  the  streamwise  variable  becomes  'tlme-llke*  and  we  are  able  to  compute  the 
invlscid  field  from  the  equations  presented  here.  This  approach  readily  admits  the  treat¬ 
ment  of  norvaxisymmetrlc,  for  example  square  cross-section,  bodies. 
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Various  shapes  which  result  from  the  Newton-Busemann  formula  shall  be  discussed  in 
the  next  section. 

As  mentioned  earlier,  for  the  purposes  of  this  research,  the  shapes  so  derived  were 
assumed  to  be  associated  with  the  solid  body  plus  its  viscous  displacement.  The  validity 
of  the  assumption  in  hypersonic  flow  rests  on  the  fact  that  the  displacement  thickness, 
6*,  approximately  equals  the  boundary  layer  thickness,  6,  and,  from  considerations  of 
continuity,  the  slope  of  the  inviscid  streamlines  at  the  edge  of  the  boundary  layer  equals 
the  streamwlse  gradient  of  the  displacement  thickness: 

tan  ©  s  v  /u_  =  d6*/dx  . 

Therefore,  in  hypersonic  flow,  the  edge  of  the  boundary  layer  is  approximately  an  invi¬ 
scid  flow  streamline.  It  Is  well  known  that,  in  inviscid  flow,  any  streamline  may  be 
replaced  by  a  solid  boundary.  As  a  result,  the  solid  body  shapes  which  support  the  com. 
puted  optimum  viscous  displacement  shapes  may  be  obtained  simply  by  subtracting  6*. 

Several  caveats  apply  here,  however.  First,  if  the  flowfieid  conditions  do  not  justify  the 
approximation  that  6*  <*  4*  the  6*  correction  procedure  must  be  modified.  In  addition,  we 
note  that  in  hypersonic  flow  this  approach  is  not  equivalent  to  the  heuristic  use  of  6* 
corrections  in,  for  example,  internal  flows  where  the  equivalent  uniform  inviscid  flow  is 
often  introduced  on  the  basis  of  global  conservation  of  mass  flow  arguments.  Finally, 
the  definition  of  6*  for  slender  bodies  must  be  treated  with  some  care  as  explained  In 
the  following  article. 
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D  3  local  body  diameter 

H  =$hapetactor,6*/$ 

k  =emnalpy 

M  -  Mach  number 

5  *  exponent  i  n  power  law  property  profiles 

Pr  =  Prandtl  number 

R  s  local  body  radius 

Re  =  Reynolds  number 

T  =  temperature 

u  s  stream  wise  velocity  component 

W  3  two-dimensional  width,  Eq.  (3) 

x  3  stream  wise  coordinate 

y  3  coordinate  normal  to  body  surface 
a  3  local  slope  of  body  surface,  Eq.  (?) 

1  *  ratio  of  specific  heats 

6  3  boundary-layer  thickness 

4*  >  displacement  thickness,  Eq.  (1) 

6*  3  enthalpy  thickness 

e  3  function  of  Prandtl  number  and  flowfield 

parameters 
0  *  density 

9  *  momentum  thickness,  Eq.  (2) 

r  *  transverse  curvature  parameter,  Eq.  ( 10) 

a  *  normalized  y  coordinate,  y/6 

—  » indicates  quantities  nondimensionalized  by  8  tot 

integral  thicknesses  and  by  frcestream  values  for 
dependent  variables 
Subscripts 

aw  m  adiabatic  wall  conditions 

o»  ■freestream  conditions 

o  *  integral  parameters  in  the  form  of  Eqs.  (1 1)  and  (12) 
2D  •  two  dimensional 

w  ■  wall  conditions 

BOUNDARY-layer  displacement  and  momentum 
thicknesses  may  be  introduced  by  conceptually 
uniforaixing  the  flowfield  properties  throughout  the 
boundary  layer.  In  order  that  the  mass  flux  past  the  wall  in  the 
streamwlse  direction  equal  that  in  the  boundary  layer  with 
nonuniform  flow,  the  wall  (and  therefore  the  external  flow 
streamlines)  must  be  shifted  a  distance  8*  normal  to  its  own 
plane.  Expressed  in  mathematical  terms,  this  can  be  written 

Similarly,  the  net  momentum  defect  of  the  nonuniform  flow 
in  the  boundary  layer  rauu  be  matched  by  the  unifeemized 
flow,  forcing  the  total  momentum  In  a  layer  of  thickness  f 
adjacent  to  the  wall  to  be  accounted  for  through  the 
relationship 

( 3  j  (2) 

In  both  Eqs,  (I)  and  (2)  a  coordinate  system  which  is  at¬ 
tached  to  the  wall  has  been  utilized  and  the  y-axis  is  assumed 
to  be  normal  to  the  local  body  surface.  In  the  two- 
dimensional  configuration,  the  differential  area  through  the 
boundary  layer  is 

(!) 

where  Ik' is  f<he  local  transverse  flow  field  dimension.  With  this 
expression  applied  to  Eqs.  (t)  and  (2).  the  following  classical 
formulas  can  be  derived 

**W-j*(/“*d>de  (4) 

and 

d^wj^/^uide  01 


The  terms  on  the  right  hand  side  of  Eqs.  (4)  and  (S)  are  also 
the  exact  form  of  the  grouping  of  variables  which  appear  in 
the  two-dimensional  integral  boundary-layer  equations. 

For  the  axisymmetric  situation,  the  local  differential  area 
through  the  boundary  layer  may  be  written 

dA=2r(R  ±y  cos  a)dy  (6) 


where: 


ctwtan''{d/7/dx)  (7) 

The  +  or  -  sign  in  Eq.  (6)  is  applicable  to  external  or  internal 
flows  respectively  and  R  is  the  local  perpendicular  distance 
from  the  axis  of  symmetry  to  the  body  surface.  Utilizing  this 
differential  area  in  Eqs.  (!)  and  (2)  results  in  the  following 
expressions  for  the  displacement  and  momentum  thickness 

TF*+ST-{^U-*m/+2T4)d4»0  (8) 


and: 


#+8-J^W(i-i))(/+2r*)dt»0 


(9) 


Two-dimensional  flows  are  embedded  in  this  general  for¬ 
mulation  as  the  limiting  case  of  r«0,  where 

xm  ±8  cos  a  ID  (10) 


and  the  a  sign  carries  the  same  significance  as  previously 
noted.  Since  r,  4*  and  i  are  significantly  less  than  unity  for  the 
majority  of  cases  of  Interest,  Eqs.  (8)  and  (9)  can  generally  be 
approximated  by 

4*>j^/-pu)(/+2r,)dg  (11) 


and 


I,  »{/»</-«)(/+*•) d«  (12) 


Upon  integrating  the  axisymmettk  equations  of  motion 
across  the  boundary  layer  (Ref.  I),  the  resulting  momentum 
integral  equation  explicitly  involves  the  terms  and  4*lo  the 
torm 


111) 


where 


Hn8\t9<  (14) 

and  4;  and  t*  (i.e.  and  9  J)  are  calculated  from  the 
dimensional  counterparts  of  Eqs.  ( 1 1)  and  (12). 

it  it  now  apparent  that  the  uniformly  valid  physical 
descriptions  of  the  displacement  and  momentum  thickness,  as 
expressed  mathematically  by  Eqs.  (I)  and  (2).  supply 
relationships  which  are  identical  to  those  appearing  in  the 
governing  momentum  integral  equations  for  values  of  r 
generally  encountered.  However,  for  cases  of  large  r  with 
flowfield  conditions  which  generate  targe  values  of  4*  and  9, 
Eqs.  (I)  and  (9)  should  be  applied  when  calculating  the  in¬ 
tegral  thickness  parameters  which  have  been  described 
physically  (4*  9),  while  the  appropriate  forms  of  4*  and  9„  are 
utilized  in  the  solution  of  the  equations  of  motion. 

The  principal  point  of  this  Note  js  to  emphasize  those 
situations  in  which  the  terms  8*  and  9  are  being  considered, 
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apart  from  their  relationship  to  terms  appearing  in  the 
governing  conservation  equations,  and  in  which  the  ex¬ 
pressions  generally  employed  for  their  evaluation  are  those 
appropriate  to  5*0  and  8„.  For  example,  consistent  empirical 
correlation  of  the  skin  friction  coefficient  with  momentum 
thickness  and  the  shape  factor  with  the  various  independent 
flowfield  parameters  should  extend  more  universally  into  the 
large  r  domain  utilizing  the  general  expressions  contained  in 
Eqs.  (8)  and  (9).  In  addition,  the  use  of  displacement  thickness 
distributions  to  generate  “effective  bodies”  in  hypersonic 
flow  utilizing  the  tangent  cone  approximation  or  in  the  study 
of  the  entry  region  of  an  axisymmetric  pipe  flow  would 
benefit  from  the  use  of  Eqs.  (8)  and  (9)  which  result  directly 
from  the  simple  physical  descriptions  of  displacement  and 
momentum  thicknesses. 

The  results  of  the  investigation  reported  here  provide 
general  indications  of  the  limits  on  r  for  which  the  terms  8a, 
P„  and  //,,  that  arise  from  the  integral  conservation 
equations  are  valid  measures  of  8,  P  and  H  which  have  their 
basis  in  the  development  leading  to  Eqs.  (1 )  and  (2).  Given  the 
dependence  of  8,  l *  and  H  on  r  and  the  parameters  which 
characterize  the  flowfield,  the  values  of  80,  and  H0  can  be 


obtained  from  the  following  relations,  derivable  from  Eqs. 
(8).  (9),  (11).  and  (12): 

(15a) 

p 

(15b) 

H  /<*tT  ' 

USc) 

These  equations  also  provide  the  percentage  difference 
between  the  two  sets  of  integral  thickness  variables  of  interest. 
In  order  to  detcrmtnc  the  quantitative  effect  of  r  on  the 
displacement  and  momentum  thicknesses  for  a  representative 
case  of  tmerest,  the  following  formulas  for  the  distributions 
of  velocity  and  temperature  through  a  constant  pressure 
turbulent  boundary  layer  have  been  utilized: 

D»,'  ' 

(16a) 

tvo  +  bl'  v  +  tV'* 

(16b) 

For  these  computations,  the  term » in  the  coefficient  6  of  Eq. 
(16b).  which  characterizes  the  lack  of  similarity  between  the 
velocity  and  temperature  fields  for  Praodtl  numbers  different 
from  one.  has  been  set  equal  to  unity. 

If  it  is  assumed  that  the  pressure  is  constant  across  the 
boundary  layer  (another  aspect  of  targe  r  layers  to  be 
quantitatively  evaluated),  Eqs.  (16)  can  be  applied  to  Eqs  (t) 
and  (9)  and  the  resultant  integral  parameters  computed  as  a 
function  of  the  principal  independent  flowfield  variables 

*».«(  >,  Pt,  r.,.W»V,t )  (17) 

The  effects  of  i.Pt.M*  and  A’ did  not  significantly  affect  the 
specification  of  a  characteristic  r  for  which  she  integral 
boundary-layer  characteristics,  as  evaluated  b>  Eqs.  (6)  and 
(9).  depart  from  the  value)  calculated  from  Eqs.  o  1  j  and  (12). 
Therefore,  ai  M,  *  2  and  A'«  7,  the  influence  of  t .  and  *  on 
the  integral  parameters  f  and  H  hat  been  computed  for  air 
with  >*1.4  and  fv*.  72  as  shown  in  Fig.  I  The  data  en¬ 
compass  a  w  ide  range  of  values  for  f  and  T « ,  net  because  they 
are  physically  likely  to  occur  but  solely  to  expose  the  variable 
dependences  of  interest. 

Mach  numbers  from  0-4  and  power-law  exponents  from 
J- 10  at)  eshibiied  curses  of  similar  shape  As  higher  Mach 
numbers  the  sensitivity  of  H  to  r  for  cold  wail  conditions 


Fig.  1  Dependence  of  t  tad 
H  on  t  and  Tw  for  t  constant 
pressure  turbulent  boundary 
layer.  =  2,  A'«1,  >*1.4 
and  Pr*>  0.72. 


becomes  increasingly  positive  white  the  shift  in  position  of  the 
curves  for  a  varying  N  amounted  to  a  few  percent  for  each 
integer  change  in  the  exponent. 

Utilizing  the  results  presented  in  Fig.  1  and  Eqs.  (15).  the 
percentage  difference  between  the  two  displacement 
thicknesses  at  r*  1.0  and  t.  -  4.0  has  been  calculated  to  be 
40T*  while  the  shape  factors  differ  by  28<T*.  It  is  thus  dear 
that,  for  high  r  boundary  layers,  utilization  of  the  physically 
describablc  entities  P  and  8  must  be  carefully  distinguished 
from  the  terms  Pu  and  8?  which  arise  in  the  governing  in¬ 
tegral  conservation  equations. 

To  this  point,  other  commonly  used  integral  boundary- 
layer  thicknesses  (Ref.  3)  have  been  neglected,  inasmuch  as 
their  usefulness  lies  primarily  in  the  fact  that  they  allow  the 
governing  equations  to  be  presented  in  a  compact  for¬ 
mulation.  Their  application  to  procedures  outside  the  scope  of 
solving  the  equations  is  relatively  limited  One  infrequent 
exception  is  the  use  of  the  enthalpy  thickness  (6* )  to  correlate 
the  dependence  of  the  Stanton  number  on  the  Reynolds 
number.  Re**,  as  shown  In  Ref.  2. 

The  enthalpy  thickness  can  be  introduced  on  a  purely 
physical  basis  in  a  manner  identical  to  that  used  for  the 
momentum  thickness.  The  resulting  mathematical  expression 
for  this  physical  entity  is 


M+K-  ]  w<tf-/)(/+2rq)de*0  (It} 

By  integrating  the  axisymmetric  energy  equation  through  the 
boundary  layer,  the  term  arises  naturally  from  the  in¬ 
tegration  procedure,  where 

w  (W~m/+2r,)de  (19) 

w 

Eq.  (I  I  Jean  therefore  be  written: 

(20) 

However,  for  a  perfect  gas  boundary-lam  flow  wiih  constant 
pressure  across  the  sttcous  region,  it  is  possible  »e  show  that 

**'“**“L  (21) 


If  a  velocity  profile  or  the  form  presented  in  Eqs.  (16)  is 
utilized.  Eq.  (21)  results  in 


[A*  t/  +  rn.v+/n 
l  tA*/){2N+/f  J 


(22) 


Therefore,  the  generalised  enthalpy  thickness,  with  which  the 
convective  heat  transfer  is  to  be  correlated,  may  be  computed 
from  Eq.  (20).  written  in  the  form 


v  /  Verf)  /it  (2» 


600 


AIAA  JOURNAL 


VOL.  15,  NO.  4 


with  hho  being  evaluated  utilizing  Eq.  (22)  and  the  results  of 
computations  similar  to  those  presented  in  Fig.  1 . 

Finally,  the  generalized  definitions  for  each  of  the  integral 
thicknesses  employed  here  are  compatible  with  the  limiting 
case  of  axial  boundary-layer  flow  over  a  needle.  In  this 
situation  r  approaches  infinity  and  each  of  the  integrals 
reduces  to  the  form  of  the  following  equation  for 
displacement  thickness 

5i”>i  =  (S!a“^)d’’2)  *  <24> 
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With  this  justification,  the  solid  bodies  which  support  the  optimum  viscous  displacement 
bodies  may  be  established  by  first  computing  the  displacement  thickness  and  then  sub¬ 
tracting  it  from  the  optimum  snape.  This  procedure  is  necessarily  iterative  inasmuch  as, 
although  the  inviscid  pressure  field  is  known,  the  solid  surface  is  not.  Therefore,  the 
location  of  the  surface  boundary  conditions  is  also  unknown.  The  iteration  begins  by 
assuming,  initially,  that  the  solid  body  shape  is  the  computed  optimum  shape.  The  resul¬ 
tant  boundary  layer  displacement  is  then  subtracted  from  the  optimum  shape  and  the 
entire  procedure  is  repeatec.  At  high  Reynolds  number,  convergence  is  obtained  within 
a  few  cycles.  At  low  Reynolds  number,  convergence  may  be  ensured  only  by  partially 
correcting  the  first  guess  solid  body  sh  .pe  so  as  to  ensure  that  a  large  overshoot  does 
not  occur. 

Computation  of  the  viscous  flow  displacement  and  skin  friction  at  each  Iteration  Is 
carried  out  within  the  framework  of  the  following  equations,  utilizing  standard  notation. 
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AX  I  SYMMETRIC  GOVERNING  EQUATIONS  OF  NOTION 
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Cnee  the  solid  body  shapes  have  been  computed,  the  internal  volume  distribution  may 
be  calculated.  During  Phase  I,  a  simple  right  circular  cylindrical  penetrator  shape  was 
assumed.  The  maximum  volume  of  such  a  cylinder  which  can  fit  within  the  solid  body 
shape  divided  by  the  total  body  volume  is  defined  as  the  prismatic  coefficient,  P.  During 
subsequent  work  other  reference  penetrator  shapes,  including  dual  in-line  cylinders  of 
different  diameter,  may  be  considered. 

The  loss  in  kinetic  energy  along  the  trajectory  may  then  be  computed  from  simple  rigid 
body  dynamics: 


V/V0  =  exp  (-0  XC^/2) 

where:  V0  =  launch  velocity 

Cfj  ■=  a^eiage  along  x 

x  =  x/:_r|;f 

X  =  length  of  flight 

'-ref  =  total  volume/base  area 

p  =  ambient  density/average  projectile  density 

For  preliminary  computations,  the  average  projectile  density  may  be  approximated  as 
the  prismatic  coefficient,  P.  times  the  density  of  the  penetrator.  Therefore,  ('or  a 
given  length  of  flight, ambient  density  and  penetrator  material  density,  the  minimum 
velocity  degradation  Is  associated  with  minimizing: 

co^ 
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It  is  clear  that  Lfe#  varies  from  shape  to  shape  so  that  the  actual  degradation  in  velo¬ 
city,  even  for  a  given  range  and  equal  CD/P,  changes  also. 

Unfortunately,  the  competition  between  candidate  shapes  is  not  quite  so  straightforward 
when  real-life  operational  considerations  are  introduced.  For  example,  the  launch  device 
'barrel'  may  constrain  the  base  diameters  to  be  equal.  In  addition,  we  may  insist  on  com. 
paring  projectiles  with  the  same  total  launch  mass.  Again  neglecting  the  residual  projec¬ 
tile  structure,  the  total  mass  may  be  written: 

m  =  A  L  f  P 

pen.  ref 

If  only  the  base  areas  of  candidate  projectiles  must  be  equal,  their  masses  will  vary  as 
the  product  of 


for  the  same  penetrator  density.  If,  however,  the  masses  must  also  be  equal,  then  the 
product: 


must  be  constant.  This  may  be  effected  by  varying  the  aspect  ratio  of  the  computed 
optimal  configurations  from  shape  to  shape*.  Clearly,  aspect  ratio,  mass  and  base  area 


•In  so  doing,  however,  projectile  stability  considerations  for  a  given  spin  or  fin  geometry 
also  restrict  the  acceptable  range  of  aspect  ratio. 
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may  not  all  be  constrained.  (A  somewhat  artificial  way  out  of  this  problem  is  to  assume 
that  any  ullage  will  be  filled  with  exactly  that  amount  of  (non-penetrator)  mass  required 
to  make  the  payload  masses  equal.)  For  purposes  of  the  Phase  I  work,  it  is  sufficient  to 
require  the  masses  and  aspect  ratios  to  be  equal  for  all  candidate  shapes.  The  base  areas 
shall  then  vary  inversely  as  the  product: 


As  a  result,  any  launch  device  'barrel'  diameter  constraint  Is  Ignored. 

The  final  'real-life'  problem  is  the  requirement  to  compute  Cp  at  off-design  Mach  and 
Reynolds  numbers  in  order  to  determine  its  average  value  along  the  trajectory. 

Given  the  ability  to  predict  the  momentum  and  kinetic  energy  of  the  projectile  at  the 
target  location,  the  details  of  the  Initial  stages  of  impact  are  computed  from  the 
Navier-Stokes  equations  utilizing  physical  properties  appropriate  to  the  materials  in- 
volved*'*: 
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TRANSIENT  MULTI-DIMENSIONAL  EQUATIONS  OF  MOTION 
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Following  the  work  of  MacCormack\  the  momentum  transfer  to  the  target  may  be 
computed  in  order  to  quantify  the  target  damage.  Typical  results  for  a  cylindrical 
penetrator  are  discussed  in  the  next  section. 


-16- 


Results  and  Conclusions 


On  the  basis  of  the  Newton-Busemann  pressure  equation,  the  optimal  viscous  displace¬ 
ment  body  configurations  for  five  different  cases  were  obtained.  The  results  are  presen¬ 
ted  in  the  following  figures  for  an  overall  aspect  ratio  of  five.  Since  the  projectile  mass 
is  constant  from  case  to  case,  the  base  area  is  assumed  to  vary  inversely  as  P  x  Lref 
for  a  given  penetrator  shape  and  density. 

Once  the  displacement  body  shape  was  obtained,  state-of-the-art  turbulent  boundary 
layer  computations  were  carried  out  in  an  iterative  mode  for  the  following  conditions: 

M  =  5.0  , 

Re^  *  106  ,  and 
Smooth  adiabatic  wall. 

As  a  result  of  the  Iteration,  a  consistent  combination  of  6*  and  body  radius  whose  sum 
equalled  the  computed  optimum  shape,  was  converged  upon.  Due  to  time  and  cost  cons¬ 
traints,  lower  Reynolds  numbers  and  thicker  boundary  layers  were  not  considered  due 
to  their  slow  rate  of  convergence. 

The  resultant  CD  (inviscld  plus  viscous  on  the  forebody)  and  P  are  shown  on  each  figure 
with  which  the  velocity  degragation  may  be  assessed.  In  order  to  account  for  fin  stabili¬ 
zation,  an  additional  CD  increment  is  required.  For  spin  stabilization,  the  aforementioned 
boundary  layer  computation  has  been  modified  to  allow  the  displacement  thickness  sub¬ 
traction  process  and  skin  friction  computation  to  be  carried  out.  Typical  results  for  the 
conical  body  at  two  spin  rates  ate  presented  in  Figures  7a  -  7d. 
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For  an  aluminum  penetrator,  with  a  length  to  diameter  ratio  of  one,  impacting  a  semi¬ 
infinite  aluminum  target,  typical  results  of  a  Navier-Stokes  simulation  of  the  initial 
stages  of  impact  are  depicted  in  Figure  8.  On  the  basis  of  the  various  fluxes  of  momenta, 
target  damage  effectiveness  may  be  quantified. 

As  a  result  of  the  Phase  I  work,  it  is  clear  that  a  systematic  consideration  of  (inviscid 
and  viscous)  drag,  internal  volume  distribution  and  impact  effectiveness  is  feasible  in 
order  to  establish  optimal  hypervelocity  projectile  shapes.  Subsequent  wind  tunnel  valida¬ 
tion  of  the  results  from  an  extended  application  of  this  analysis  is  required  orlor  to  the 
execution  of  ballistic  range  and  impact  experiments. 
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Table  1. 


Minimum  Drag  Projectile  Design  Cases  Studied 


Case  1.:  Figures  2a  &  2b:  Given  length  and  base  diameter. 


Case  2.:  Figures  3a  &  3b:  Given  length  and  volume. 


Case  3.:  Figures  4a  &  4b:  Given  length  and  surface  area. 


Case  4.:  Figures  5a  &  5b:  Given  base  diameter  and  volume. 


Case  5.:  Figures  6a  &  6b:  Given  base  diameter  and  surface  area. 


Case  6.:  Figures  7a  -  7b:  Case  6  with  two  spin  rates 
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Figure  1. 


Schematic  of  Optimization  Considerations 
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Case  1. 


Minimum  Drag  Shape  for  a  Given  Length  and  Base  Diameter 

Figure  2a.:  Original  and  corrected  shapes. 

Figure  2b.:  Skin  friction  and  pressure  distributions. 
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Case  2. 


Minimum  Drag  Shape  for  a  Given  Length  and  Volume 

Figure  3a.:  Original  and  corrected  shapes. 

Figure  3b.:  Skin  friction  and  pressure  distributions. 
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Fig.  3a. 
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Minimum  Drag  Shape  for  a  Given  Length  and  Surface  Area 


Figure  4a.:  Original  and  corrected  shapes. 

Figure  4b.:  Skin  friction  and  pressure  distributions. 
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Fig.  4a. 
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Case  4. 


Minimum  Drag  Shape  for  a  Given  Base  Diameter  and  Volume 


Figure  5a.:  Original  and  corrected  shapes. 


Figure  5b.:  Skin  friction  and  pressure  distributions. 
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Fig.  5a. 
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Fig.  5b. 
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Case  5. 


Minimum  Drag  Shape  for  a  Given  Base  Diameter  and  Surface  Area 


Figure  6a.:  Original  and  corrected  Shapes. 

Figure  6b.:  Skin  friction  and  pressure  distributions. 
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Fig.  6a. 
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Case  6. 


Case  5.  Repeated  with  Body  Rotation 


Figure  7a.:  Original  and  corrected  shapes  at  a  spin  rate  of 
2  x  10s  rpm. 


Figure  7b.:  Skin  friction  and  pressure  distributions  at  a  spin 
rate  of  2  x  io5  rpm. 


Figure  7c.:  Original  and  corrected  shapes  at  a  spin  rate  of 
1  x  10s  rpm. 


Figure  7d.:  Skin  friction  and  pressure  distributions  at  a  spin 
rate  of  1  x  io5  rpm. 
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Fig.  7a. 
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Fig.  7c. 
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Figure  8.  Navier-Stakes  Impact  Solution  For  Cylindrical  Penetrator 
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Unified  Inlet/Diffuser  Design  by  an  Inverse  Method 

G.  J.  Hokenson* 
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The  equations  governing  the  two-dimensional,  steady,  inviscid  flow  of  an  incompressible  fluid  have  been  for¬ 
mulated  in  an  inverse  manner  to  obtain  the  resultant  shape  of  semi-flush  inlet/curved  diffusers  when  the  surface 
velocity  distribution  Is  specified.  The  results  of  calculations  illustrating  the  effect  of  each  of  seven  principal 
design  variables  on  a  nominal  inlet  configuration  are  presented.  One  of  the  configurations  obtained  was  chosen 
for  experimental  testing  to  compare  the  on-  and  off-design  performance  of  inlets  generated  by  this  procedure 
with  that  of  inlets  obtained  by  more  conventional  methods.  General  results  and  applications  of  this  method  are 
discussed. 


Nomenclature 

AR  =area  ratio,  h,/.h< 

d  = vertical  distance  from  inlet  lip  to  ramp  surface  at 
-«MFig.  1) 

DF  =drop  fraction,  d/h, 

h,  =  vertical  height  of  ramp  above  inlet  lip.  (Fig.  1 ) 

h  =  inlet  height  (Fig.  I) 

IVR  =inlet  velocity  ratio,  U(/Um 

K  =  inlet  loss  coefficient,  (Pr  ~PT  )/q  where: 

q^qfot  IVR>  1,4=0,  for/kfecl. 

L  =  ramp  length  (Fig.  I) 

Ld  = diffuser  length  (Fig,  l) 

LVR  -  lip  velocity  ratio,  (/,,  /  Um 

In  -  natural  base  logarithm 

Pr  ~  stagnation  pressure 

q  adynamic  pressure 

SM  ajtreamwise  position,  measured  in  inlet  heights  up¬ 
stream  from  the  lip.  at  which  the  ingested  stream- 
tube  departs  from  a  line  which  parallels  the  ramp 
(Fig.  1) 

u  »  component  of  velocity  in  *  direction 
U  ■  velocity  vector  magnitude,  (u-  +V )  ti 

y  «  component  of  velocity  in  y  direction 

x  a  horizontal  Cartesian  coordinate  normalized  by  L 
with  origin  at  end  of  upstream  uniform  flow  region 
y  » vertical  Cartesian  coordinate  normalized  by  L  with 

arbitraryorlgin 

a  •  velocity  vector  angle,  tan '' (t'/u) 

v  it)  *  Laplaclan  operator  in  «  -  coordin'  ‘es 

4  "b’/d, 

*iO  a  finite  difference  grid  spacing  in  o 
A*  a  finite  difference  grid  spacing  in  <!> 

0  a  velocity  potential:  Eq.  (5)  and  (6) 

o*  «o-o, 

vi  a  stream  function;  Eq.  (5)  and  (6) 

r  a  ingested  streamtube  length  to  height  ratio,  L(hm, 
(Fig.  I) 

r0  « diffuser  length  to  height  ratio.  Ln.'h, 

$  *  lip  angle  with  respect  to  the  vertical  (Fig.  I ) 

w  •'h.tli* 
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Subscripts 

e  = conditions  at  diffuser  exit 
/  =  streamline  index 

o  =conditionsato(Fig.  2) 

r  =conditions  at  r  (Fig.  2) 

i  =conditionsati(Fig.  2) 

»  =  conditions  along  ramp  surface  infinitely  far  up¬ 

stream 

t  -  conditions  at  inlet  lip  (Fig.  3) 

/  «  lowest  bounding  streamline  in  computation  region 

(Fig.  2) 

2,3  =dividing  streamline  which  splits  at  the  lip  (Fig.  2) 

4  *  streamline  which  forms  ramp  and  upper  surface  of 

diffuser  (Fig.  2) 

Introduction 

NTAKE  ducting  in  which  an  opening  along  a  solid 
bounding  surface  io  the  flow  field  is  used  to'draw  off  fluid 
has  been  put  to  so  many  applications  as  to  prohibit  their  ex¬ 
plicit  enumeration.  Of  prime  importance  among  these  ap¬ 
plications  are  engine  air  inlets,  auxiliary  equipment  inlets  and 
boundary -layer  suction  slots,  whether  used  in  wind  and  water 
tunnels  to  "cleanse"  the  flow  or  to  retard  separation  in  a 
variety  of  adverse  pressure  gradient  situations.  For  the 
majority  of  these  inlet  applications,  a  coupled  diffuser  is  in 
evidence  which  allows  the  system  to  operate  efficiently  with¬ 
in  the  power  constraints  of  the  device  which  sustains  the 
flow.  Possibly  equally  as  common  is  to  have  the  diffuser  curve 
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inboard  from  the  flow  surface  along  which  the  fluid  is  being 
drawn  off.  Initial  experimental  results  obtained  by  the 
authors  on  inlets  such  as  those  shown  in  Fig.  1  indicate  that 
the  inlet  geometry  and  flow  conditions  upstream  of  the  lip, 
and  to  a  finite  but  lesser  extent  on  the  downstream  external 
lip,  are  strongly  coupled  to  the  pressure  recovery  and  overall 
character  of  the  flow  inside  the  curved  diffuser.  This  is  not  a 
particularly  surprising  fact,  yet  it  is  one  which  has  received  lit¬ 
tle  intensive  study.  In  general,  the  interaetion  effect  of 
series/ parallel  hydrodynamically  coupled  components  in  a 
fluid  system  on  the  detailed  flowfield  is  not  yet  understood.  In 
engineering  applications,  established  practice,  guided  by  ex¬ 
perience  and  some  excellent  classical  experimental  studies 
(e.g.,  Ref.  4-6),  initiates  the  general  design.  Utilizing  this 
procedure,  many  very  efficient  inlets  have  been  developed  af¬ 
ter  considerable  empirical  iteration  on  the  design  has  been 
completed  and  which  essentially  compensates  for  these 
unknow  n  interactions. 

Dominant  in  the  open  literature  with  respect  to  their  scope 
and  intensity  of  experimental  results  are  the  publications  of 
Mossman  (e.g.,  Ref.  4)  which,  for  inlet  geometries  of  the 
general  type  investigated  here,  describe  an  optimum  external 
pre-diffusion  (indicated  by  an  1VR  less  than  one)  for  a  given 
set  of  variations  of  a  specific  inlet  configuration.  These  results 
have  proved  to  be  an  invaluable  guide  in  the  development  of 
inlets  with  relatively  simple  attached  diffusers;  however,  the 
introduction  of  more  complicated  large-area-ratio,  curved 
diffusers  changes  this  picture  significantly.  The  generation  of 
curved  diffuser  shapes  with  relatively  large  expansions  has 
been  well  documented  in  Ref.  3,  supported  by  the  analytical 
background  developed  in  Ref,  I .  The  excellent  experimental 
data  obtained  at  Stanford' 3  has  served  as  a  guide  to  the  ex¬ 
tent  to  which  the  boundary  layers  on  each  wall  of  a  curved  dif¬ 
fuser  can  support  adverse  pressure  gradients,  and  indirectly 
the  range  of  geometrical  parameters  which  will  result  in  an  ef¬ 
ficient  pressure  recovery  device.  There  is  little  information 
available,  however,  which  can  serve  as  a  guide  to  the  unified 
design  of  an  inlet/diffuser  combination  when  large  amounts 
of  pre-diffusion,  varying  angle  of  incidence  of  the  dividing 
streamline  onto  the  inlet  lip.  and  external  lip  geometry  all  in¬ 
teract  with  the  internal  diffuser  area  ratio,  curvature,  and 
length-io-height  ratio.  This  gap  tn  the  information  available 
in  the  open  literature,  coupled  with  the  aforementioned  ex¬ 
perimental  results,  has  directed  the  authors  to  look  at  some 
simple  analytical  tools  for  use  in  the  design  of  such  inlet 
systems  and  compare  them  experimentally  to  configurations 
which  were  designed  without  benefit  of  the  analysis  presented 
here. 

With  regard  to  the  direct  analysis  of  inlets  and  diffusers 
with  a  specified  geometric  configuration,  .the  current  state-of- 
the-art  allows  the  full  equations  of  motion  with  a 
sophisticated  Reynolds  stress  tensor  model  to  be  applied  to 
the  duet  This  approach  is  not.  however,  an  efficient  wall  con¬ 
tour  design  technique  and  generally  should  be  employed  to  ex¬ 
pose  the  nonideal  fluid  flow  phenomena  associated  with  a 
family  of  optimum  inlet  shapes  generated  by  an  inverse 
method  such  as  the  one  employed  here.  To  the  authors’ 
knowledge,  no  such  calculations  utilizing  the  full  equations  of 
motion  have  been  conducted  on  the  generic  type  of  tnlet*  dif¬ 
fusers  considered  here  However,  some  recent  experimental 
studies  carried  out  at  NASA  Lewis'  evaluated  annular  Grif¬ 
fith  diffusers  whose  design  was  guided  by  the  results  of  Ref.  8. 
In  this  design  procedure,  the  diffuser  wall  contour  is  obtained 
analytically  whereas  the  details  of  the  boundary-layer  suction 
slot  region  arc  not  emphasized  The  results  ot  Ref.  7  indicate  a 
significant  itnproscmcm  m  the  performance  of  short  annular 
diffusers  with  high  suctum  flow  rates.  On  the  basis  of  the 
studies  which  are  reported  here,  the  detailed  contours  of  the 
boundary -layer  cutolt  slot  and  associated  ducting  could 
possibly  be  optimized  to  improve  diffuser  recovery  at  a 
snitiller  suction  How  rate.  It<  support  of  this  hypothesis,  it  has 
been  the  authors*  experience  that,  while  observing  a  hot  *yt 


anemometer  signal  in  a  boundary  layer  downstream  of  a  suc¬ 
tion  slot,  the  entire  character  of  the  flow  is  extremely  sensitive 
to  the  inlet  IVR ,  lip  position,  and  lip  angle  for  a  given  overall 
suction  flow  rate. 


Analytical  Development 

The  design  of  a  semiflush  (i.e.,  DF> 0)  inlet/curved  dif¬ 
fuser  has  been  approached  in  an  inverse  manner  to  obtain  the 
resulting  physical  shape  of  the  bounding  surfaces  to  the 
flowfield  when  the  surface  pressure  or  velocity  distribution  is 
specified.  The  equations  governing  the  two-dimensional,  in¬ 
compressible,  inviscid  flow  appropriate  to  this  model  of  the 
flowfield  have  been  solved  over  the  entire  region  associated 
with  the  inlet,  thus  resulting  in  a  unified  shape  which  is  totally 
consistent  with  any  local  change  in  the  specified  surface  con¬ 
ditions.  The  calculation  region  includes  the  inlet  ramp, 
curved  diffuser,  external  lip,  and  the  remote  streamline  which 
is  the  lower  surface  of  the  computation  region  shown  in  Fig. 
2. 

The  governing  equations  of  continuity  and  irrotationality 
expressed  in  spatial  Cartesian  coordinates  and  dimensional 
physical  variables  are 


Continuity 


and 

irrotationality 


bu  dv 
dx  by 


=  0 


bu  b  v 
by  dx 


(1) 


(2) 


The  dependent  variables  can  be  transformed  from  the  x  and  y 
components  of  velocity  to  the  velocity  vector  magnitude  and 
angle  via  the  following  definitions 

U~\u'  +t>-  (3) 


and 


<z«tan"'(i»/u)  (4) 

In  addition,  the  independent  variables  can  be  transformed 
from  a  spatial  Cartesian  coordinate  system  to  an  orthogonal 
velocity  potential-stream  function  system  via  the  following 
transformations 


do  dyl 

dx  "  ~by 


(5) 
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d4>  dd 

dy  dx 


where  the  initial  coordinates  of  each  streamline  are  obtained 
by  integrating  Eqs.  (14)  and  (15)  along  the  line©  =0  from  the 
plane  of  symmetry  streamline  to  the  streamline  of  interest  as 
shown  in  the  following  equations 


Upon  applying  these  transformations  to  the  dependent  and 


independent  variables  of  Eqs.  (1)  and  (2),  the  governing 
equations  become 

*,(o)=-j 

(18) 

Continuity 

(19) 

and 


dlnU  da 
"do"  +  3d 


Irrotationality 


It  is  now  apparent  that,  to  obtain  the  coordinates  of  each 
point  on  the  streamlines  of  interest,  the  distribution  of  a 
along  those  streamlines  must  be  obtained  from  the  known 
InU  field.  From  the  chain  rule  we  can  write 


din U  3a 
3d  3© 


(8) 


By  cross  differentiating  and  adding  Eqs.  (7)  and  (8),  we  obtain 
Laplace’s  equation  for  )nU 

vj„.lnU=0  (9) 


d©  + 


da 

3d 


(20) 


which,  utilizing  the  conditions  of  continuity  and 
irrotationality  presented  in  Eqs.  (7)  and  (8),  may  be  rewritten 


da= 


dlnU 

3d 


d©- 


dlnU 

d<t>~ 


dd 


(21) 


After  imposing  the  appropriate  boundary  conditions  on  In U 
in  the  o-d  plane,  the  entire  field  may  be  solved  (in  this  paper 
using  successive  over -relaxation)  for  the  distribution  of  InU 
throughout  a  grid  network  appropriate  to  the  inlet  geometry 
as  shown  in  Fig.  2.  Once  the  distribution  of  inU  is  known 
within  the  field  of  interest,  the  physical  shape  of  the  bounding 
surfaces  may  be  obtained  from  the  inverse  transformation  to 
Cartesian  coordinates.  Because  of  the  orthogonality  of  the  ©• 
d  coordinate  system,  the  inverse  transformation  can  be 
formed  from  the  following  chain  rule  equation' 


(10) 

^•l*  l.d*+w  l.w 

(ID 

The  definitions  of  the  independent  variables 
(5)  and  (6)  allow  us  to  write 

given  in  Eqs. 

d©x>udx+t/dy 

and 

(12) 

dd“  -odjr+udy 

U3) 

By  integrating  Eq.  (21 )  along  the  particular  streamlines  d,  we 
may  determine  the  distributions  of  a  necessary  to  complete 
the  integrations  indicated  by  Eqs  (16-19).  Upon  performing 
the  indicated  a  integration  along  d,  we  obtain 

ndlnfA 

~Sif)  dq>  (22) 


Here  again,  the  initial  streamline  angles  are  obtained  by  in¬ 
tegrating  in  the  d  direction  along  ©  =0  line  from  the  plane  of 
symmetry  streamline  to  the  streamline  under  consideration. 
This  is  done  utilizing  Eq.  (2 1 )  and  can  be  written 


a,(o)- 


i’*‘  /  dint/ \ 

J  V  /  ^aU 


(10 


(23) 


This  completes  the  integrations  which  must  be  made  utilizing 
the  previously  determined  InU  field,  resulting  in  the 
specification  of  the  location  of  each  point  on  the  streamlines 
of  interest  in  physical  space. 

Inasmuch  as  the  input  to  the  calculation  is  the  surface 
velocity  or  pressure  on  the  appropriate  boundaries  of  the 
calculation  region  in  the  ©-d  plane,  appropriate  distributions 


which  can  be  solved  for  the  appropriate  derivatives  in  Eqs. 
(I0)and;il)  resulting  in 


dr 


cosa 

l  — 

U  i 


d©- 


slno 

~TT 


d  d 


(14) 


and 


iy< 


sinn 

U 


d©+ 


CQStt 

V  I* 


dd 


<«J) 


These  results  may  now  be  used  to  obtain  the  physical  coor¬ 
dinates  of  each  point  on  any  streamline  by  integrating  across 
o  along  a  constant  d, .  The  resulting  equation  for  the  or  and  y 
positions  are 

*,(©)**,(«)■*•  j  (~r)  d©  (16) 

and 

«sinn\ 

(17) 


mu 


46- 


42 


C..J.  HOKENSON  AND  F.Y.  SU 


A1AA  JOURNAL 


must  be  selected.  For  simplicity  and  ease  in  manipulating  their 
functional  form  numerically,  velocity  profiles  of  the  form 
U~  e"-'1  were  chosen  to  form  various  pieces  of  the  distribution 
on  the  respective  surfaces  of  the  computation  region  and  are 
depicted  in  Fig.  3.  Reference  3  points  out  that  one- 
dimensional  floss  considerations  indicate  that  velocity 
distributions  of  this  type  are  generated  by  straight  wall  two- 
dimensional  channels  if  a  one-segment  exponential 
distribution  is  assumed  to  exist  over  the  entire  length  of  the 
channel.  On  all  surfaces  ot  the  computation  region  not  noted 
in  Fig.  3,  uniform  floss  may  be  specified. 

The  remote  streamline  was  moved  sufficiently  far  from  the 
inlet  so  as  to  appear  mathematically  to  be  located  at  infinity. 
The  boundary  condition  applied  at  this  surface  was  that  of 
zero  norma!  derivative  which  forced  this  streamline  to  be  a 
plane  of  symmetry  between  the  inlet  flowfield  which  was 
calculated  and  a  fictitious  mirror  image  inlet  located 
equidistant  beneath  the  plane  of  symmetry  streamline.  The 
theoretical  desirability  of  having  this  plane  of  symmetry 
available  is  apparent  in  the  utilization  of  Eqs.  (18,  19,  and  23) 
when  we  attempt  to  retrieve  the  physical  shape  of  the 
bounding  surfaces.  In  a  practical  sense,  however,  the  plane  of 
symmetry  streamline  was  located  sufficiently  far  from  the 
inlet  so  as  to  have  a  negligible  effect  on  the  resulting  inlet 
shapes,  irrespective  of  the  type  of  boundary  condition  (i.e., 
uniform  flow  or  zero  normal  derivative- which  was  applied  to 
it. 

The  effect  of  smoothing  the  assumed  velocity  distribution 
has  been  investigated;  however,  no  significant  changes  in  the 
general  inlet  shapes  were  found.  Clearly,  minor  variations  in 
the  specific  wall  shapes  ensued;  however,  no  results  were 
found  which  justified  abandoning  these  straight-line  profiles 
with  the  corresponding  case  of  manipulating  the  various 
parameters  of  interest.  It  is  to  be  emphasized  that  the  choice 
of  this  functional  form  for  the  boundary  conditions  is  ar¬ 
bitrary  and  serves  only  as  a  reasonable  distribution  within 
which  the  effects  of  IVR,  diffuser  area  ratio,  diffuser  and 
ramp  length  to  streamtube  height,  diffuser  curvature,  and  ex¬ 
ternal  lip  shape  on  the  resultant  inlet  configurations  cun  be 
studied.  There  is  an  obvious  restriction  to  the  choice  of  this 
arbitrary  distribution  oi  surface  conditions.  It  arises  from  the 
transformation  itself  which  prohibits  stagnation  points  in  the 
flow  field.  Therefore,  the  calculated  shape  ol  the  inlet  lip  will 
necessarily  be  sharp.  In  actual  practice,  a  parabolic  noise  will 
be  faired  in  locally  . 

Discussion  of  the  Results 

The  analysis  which  has  been  presented  was  programmed  in 
double  precision  and  coded  in  FORTRAN  for  execution  on  an 
IBM  3702168  computer.  The  computation  grid  was  divided 
into  a  200  Ae  >  I  $0Asl  rectangular  finite-differenee  mesh 
which  was  uniform  in  each  direction,  It  was  determined  from 
numerical  experimentation  that  an  optimum  position  for  the 


lip  streamline  occurred  at  a  yS  value  which  allowed  25%  of  the 
total  number  of  i  lines  to  be  located  within  the  inlet  stream- 
tube.  With  the  lip  streamline  located  in  this  position,  the 
solution  (obtained  through  successive  over -relaxation)  relaxed 
to  its  final  state  in  the  fewest  iterations  and  was  the  most  ac¬ 
curate  (as  indicated  by  the  final  calculated  area  ratios  of  the 
inlet  from  its  uniform  upstream  to  its  uniform  downstream)  at 
the  maximum  iteration  limit  as  shown  in  Fig.  4.  Fixing  the  lip 
at  a  particular  yt  line  forced  Ayl-  to  vary  from  case  to  case  in  or¬ 
der  that  various  diffuser  configurations  might  be  calculated 
with  different  non-dimensional  design  parameters.  Utilizing 
this  uniform  rectangular  grid  whose  aspect  ratio  varied 
somewhat  with  the  particular  inlet  being  designed,  com¬ 
putational  times  showed  some  variation.  However,  a  general 
rule  was  that,  with  an  over-relaxation  parameter  of  1.9,  the 
inlet  computation  would  arrive  at  an  asymptotic  solution 
within  400  to  600  iterations  utilizing  approximately  10 
minutes  of  CW  time. 

Once  these  preliminary  setup  constraints  of  the  com¬ 
putational  procedure  were  established,  the  parameters  of  the 
inlet  surface  velocity  distributions  were  systematically  varied. 
Of  the  parameters  available,  one  was  fixed  as  being  a 
representative  value  for  a  desirable  inlet  shape.  This  was  the 
quantity  which  reflected  the  number  of  uniform  flow  grid 
space-  which  were  placed  upstream  and  downstream  of  the 
inlet.  Extensive  numerical  experimentation  established  that  25 
uniform  flow  grid  spaces  provided  an  adequate  inlet  shape  in 
that  all  calculated  surfaces  were  insensitive  to  further  uniform 
flow  additions. 

The  remaining  inlet  design  parameters  were  varied  in  a  con¬ 
sistent  manner  about  a  nominal  inlet  configuration  to  provide 
a  comprehensive  picture  of  how  they  influenced  the  resultant 
inle*  shape.  These  parameters  are:  area  ratio  [AR),  inlet 
velocity  ratio  (IVR),  external  lip  velocity  ratio  (L  VR),  diffuser 
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Fig.  7  Effect  of  L  YR  on  inlet  shape. 


Fig. «  Effect  of  IYR  on  ink!  shape  for  4R  HR  ■*  3.0. 


lower  surface  velocity  distribution  factor  (6.1,  external  lip 
velocity  distribution  factor  (4,1,  diffuser  length  to  height  ratio 
(r(>),  ramp  length  to  streamtube  height  ratio  (r),  and  external 
lip  velocity  distribution  factor  (4, )  which  was  held  at  the  con¬ 
stant  value  of  1.0  after  it  was  determined  that  LVR  and  4, 
exercised  sufficient  control  over  the  external  lip  shape,  The  ef - 
feet  of  these  parameters  on  the  surface  velocity  distribution  is 
shown  in  Figs.  I  and  3.  Apart  from  the  basic  sensitivity  of  the 
inlet  shape  to  these  input  parameters,  secondary  inlet  charac¬ 
teristics  can  also  be  obtained  from  the  calculations.  Those 
inlet  characteristics  which  are  an  output  of  the  analysis  are 
drop  fraction  ( DF ),  inlet  lip  angle  ((f),  and  the  ramp  position 
( .S'* )  at  which  the  ingested  streamtube  departs  from  a  line 
which  parallels  the  ramp,  as  shown  in  Fig.  I . 

The  resultant  inlet  configurations  which  are  generated  when 
each  of  the  seven  principal  design  variables  are  varied  about  a 
reference  state  are  presented  in  Figs.  5-10,  Although  each  of 
the  design  variables  influence  all  parts  of  the  final  inlet  shape, 
for  the  purpose  of  general  discussion,  their  principal  effects 
can  be  separated  fairly  clearly  and  general  statements  as  to  the 
characteristics  of  the  inlets  controlled  by  each  can  be  for¬ 
mulated. 

The  effect  of  4  .  which  controls  the  diffuser  lower  surface 
velocity  distribution  as  shown  in  Fig.  ),  on  the  inlet  shape  ts 
primarily  the  degree  to  which  the  diffuser  exit  curve*  mboard 
and  directs  the  flow  away  from  the  frecstream  as  shown  m 
Fig  5.  Negligible  changes  in  the  lip  angle,  drop  fraction,  and 
ramp  shape  were  observed.  However,  of  primary  importance 
are  the  relatively  flat  ramp  and  the  thin  lip  which  allows  the 
flow  to  diffuse  significantly  prim  to  being  turned;  which  are 
characteristic  of  all  inlets  presented  here. 

Control  over  the  externa!  lip  shape  Is  exerted  primarily  by 
the  parameters  6,  and  LVR  whose  influence  on  the  velocity 
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Fig.  9  Effect  of  .4  R  on  inlet  shape. 


Fig.  10  Effect  of  r  on  ink!  dupe. 


distribution  it  also  shown  in  Fig.  3.  The  resultant  shapes  for 
various  values  of  4.  and  LVR  are  shown  in  Figs.  6  and  7 
where,  apart  from  the  lip  shape  (which  can  be  changed  from 
concave  to  convex  by  a  small  change  in  either  variable),  the 
principal  effect  is  to  lower  the  diffuser  without  significant 
change  in  the  exit  angle  as  the  external  lip  flow  acceleration  is 
increased.  Changes  in  the  drop  t  raction  are  more  sensitive  to 
changes  in  LVR  than  they  arc  to  changes  in  4..  however  the 
general  characteristics  of  a  flat  ramp  in  the  neighborhood  of 
the  inlet  lip  and  a  diffuser  which  attempts  to  separate  turning 
and  diffusion  remain  intact. 

In  tsrder  to  expose  the  dependence  of  the  in1;**  :hipc  on 
IVR,  the  total  diffusion  (i.e;,  AR'IVRi  was  held  c.0nv  ini  in 
3.0.  The  results  presented  in  Fig.  8  illustra..  vW  or .ncipal 
feature*  of  these  inlets  as  IVR  and  aR  Ate  varle.'  fwxi  ef* 
feet  which  is  worthy  of  note  is  the  sensitivity  of.iipMfop  frac¬ 
tion  to  IVR,  which  was  also  ob»erscd  on  ".me  .alcuiatiwns 
conducted  for  a  fixed  diffuser  area  raiio  as  both  <ne  IVR  and 
total  diffusion  change.  The  second  principal  tea.utc  of  the 
mteis  wuh  varying  IVR  ts  the  rapid  vltnngc  so  the  charac¬ 
teristic  shape  of  the  external  lip  The  re»uhs  u>  Fig.  8  indicate 
this  quite  clearly  when  comparing  ihe  Up  shapes  tor  HR 
■  0  63  and  0J3.  Because  it  obscured  the  iiciifc.  the  lip  tor 
/*■/?«!)  '0  was  deleted.  In  the  calculation  of  mlci  xhtipex  with 
a  fixed  diffuser  area  ratio,  the  et  test  ot  H  R  on  the  external  lip 
shape  was  even  more  pronounced  than  that  shown  <n  >  us  * 

The  effect  of  changing  the  diffuser  at .  i  ratio  at  a  -.ono.tnt 
IVR  is  shown  in  Fig.  St  where  a  tairlv  complex  Jucye  ni  the 
diffuser  exit  angle  and  diffuser  upper  surface  pomimii  ate  ol*- 
*ersed  The  ramp,  external  lip.  and  dit  (user  *«-«*-;  sufiaccarc. 
however,  relatively  invariant  wuh  tR. 
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Two  design  control  variables  whose  effect  is  not  shown  ex¬ 
plicitly  on  the  \ eloctty  distributions.  in  Fig  3  are  the  diffuser 
length  to  inlet  height  ratio  (r,.,)  and  the  ramp  streamtube 
length  to  height  ratio  (r).  These  variables  enter  the  analysis  in 
the  transformation  from  the  assumed  surface  velocity 
distribution  in  physical  space  to  the  velocity  distribution  in  0 
and  allow  the  characteristic  ©  values  of  Fig.  3  to  be  deter¬ 
mined.  The  influence  of  an  increasing  tt)  on  the  inlet  shape  is 
to  a  very  good  approximation  simply  a  stretching  of  the  dif¬ 
fuser  around  the  curvature  determined  by  the  diffuser  turning 
parameter  (A. ).  Some  additional  effects  of  secondary  im¬ 
portance  were  slight  changes  in  the  lip  shape  and  the  extent  of 
the  parallel  wall,  uniform  flow  region  ai  the  dtf (user  exit. 

The  results  for  a  varying  ?  are  shown  in  Fig  10  where  the 
reference  length  used  in  the  calculations  (the  length  of  the 
ramp  to  the  point  where  ©*©.)  is  held  •  instant  and  the 
overall  inlet  decreases  in  site  as  the  height  of  the  ingested 
streamtube  is  varied.  The  itilct  contour'  presented  tn  Fig.  10 
indicate  that  both  the  lip  angle  and  the  departure  pom:  for  the 
ramp  streamtube  from  the  ramp  are  strong  functions  of  the 
parameter  r.  These  characteristics  are  even  more  apparent 
when  the  inlets  ate  te-pioued  for  a  constant  ink)  height  with  a 
varying  ramp  length. 

In  order  to  assess  the  utility  of  this  design  method  in 
producing  an  inlet  with  adequate  presvure  recovery  capability . 
a  representative  inlet  configuration  (shown  in  Fig.  I  la  and 
libs  was  selected  for  experimental  evaluation  A  small,  two- 
dimensional  indraft  wind  tunnel  was  built  for  this  purpose  A 
centrifugal  blower  was  utilised  to  draw  room  a?r  through  a 
10:1  area  ratio  inlet  contraction  which  contained  a  series 
of  filters.  How  straightening  honeycomb  and  turbulence 
damping  screens.  Inlet  diffusers  winch  were  tested  with  the 
general  configuration  *hown  in  Fix  I  were  then  inverted  into 
the  upper  surface  of  the  channel.  Through  active  control  of 
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thc  surface  roughness,  the  boundary  layer  approaching  the 
inlet  was  adjusted  to  the  desired  thickness  (approximately 
one-half  the  inlet  height  for  the  test  data  presented  here).  At¬ 
tached  to  the  exit  of  the  inlet/diffuser  was  a  rake  of  13  total 
pressure  probes  located  in  the  throat  of  a  flow  uniforming 
venturi  immediately  downstream  of  the  diffuser  exit.  Through 
utilization  of  the  total  and  static  pressure  distribution  data  in 
the  channel  downstream  of  the  diffuser  exit,  the  pressure  loss 
and  flow  rate  through  the  duct  were  monitored.  The 
magnitude  of  the  flow  rate  through  the  inlet  was  controlled  by 
means  of  valves  located  between  the  diffuser  and  a  second 
centrifugal  blower  which  was  used  to  draw  air  from  the  main¬ 
stream  flow  within  the  wind  tunnel  and  into  the  inlet. 

The  experimental  results,  presemed  in  Fig.  12,  also 
illustrate  the  overall  inlet  loss  coefficient  K  vs  HR  for  an  inlet 
designed  by  fitting  spline  curves  threugh  the  drsired  inlet  and 
exit  positions  of  the  diffuser  and  lip.  which  was  tested  in  the 
same  wind  tunnel.  The  data  are  self-evident  in  justifying  the 
ability  of  this  method  to  generate  high  total  pressure  recovery 
inlets.  Based  on  experience  with  a  number  of  simitai  inlets  in 
the  course  of  this  investigation,  it  i»  the  authors'  belief  that  lip 
angle,  drop  fraction,  and  diffuser  turning ‘expansion 
distribution  are  the  principal  factors  affecting  the  improved 
total  pressure  recovery  supplied  by  inlets  generated  by  the 
analysis  presented  here. 

Conclusion 

Inputting  tn  arbitrary  surface  velocity  or  pressure 
distribution  to  the  inlet  design  is  a  particularly  desirable  tack 
from  severs.;'  standpoints.  First,  from  the  point  of  view  of 
retarding  separation  as  long  as  possible  It  may  be  ad¬ 
vantageous  to  specify  the  surface  pressure  distribution  so  as  to 
avoid  peaks  in  the  pressure  gradient  or  to  match  different  ex¬ 
perimentally  observed  pressure  distributions  which  have  been 
determined  to  delay  separation  on  the  ramp,  diffuser  upper 
vurface,  and  diffuser  lower  surface.  Second,  for  use  in  liquids, 
the  ability  to  specify  the  surface  presvure  distribution  is  of 
considerable  assistance  in  controlling  the  cavitation  envelope 
of  the  inlet.  Third,  the  consideration  of  imposing  a  specific 
pressure  distribution  on  parts  of  the  inlet  to  arrive  at  a 
favorable  force  and  moment  balance  can  be  studied  to  deter¬ 
mine  the  extent  to  which  it  impact*  with  the  separation  and 
cavitation  performance.  In  thn  case,  the  unified  aspect  of  this 
design  procedure  allows  one  io  study  directly  the  extent  to 
which,  for  example,  changes  m  the  external  lip  lift  force 
distorts  the  diffuser  or  inlet  ramp  geometry  for  a  given  ramp 
diffuser  velocity  distribution  or  forces  a  distorted  ramp- 
diffuser  velocity  distribution  for  a  given  ranvp-dtffuser  shape. 
Fourth,  when  a  variable  geometry  inlet  must  be  designed  to 
operate  efficiently  at  a  variety  of  design  points,  the  inverse 
method  clearly  illustrates  which  surfaces  must  be  movable  and 
in  which  general  direction  they  must  move  to  allow  the  inlet  to 
be  tuned  io  its  optimum  performance  in  that  particular 
flow  field.  Filially  ,  the  extern  to  which  certain  specified  con¬ 
figurations  impose  unavoidable  How  pnenotnena  can  be 
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studied  by  this  method.  An  example  of  this  involves  the  semi- 
flush  nature  of  the  inlets  presented  here.  This  inlet  charac¬ 
teristic  is  a  consequence  of  the  fact  that  the  analysis  imposes  a 
uniform  flow  at  the  inlet  lip  by  having  both  the  diffuser  upper 
and  lower  surface  velocity  profiles  originate  at  In 
IVR  as  shown  in  Fig.  3.  It  is  equally  feasible  to  generate  flush 
inlets  by  having  the  diffuser  lower  surface  velocity 
distribution  originate  at  values  somewhat  smaller  than  In 
IVR  and  the  diffuser  upper  surface  velocity  distribution  at 
values  somewhat  larger  than  In  IVR,  thus  displacing  the  inlet 
lip  downstream  and  more  nearly  flush.  Therefore,  we  can  in¬ 
fer  that  in  a  direct  attack  upon  the  inlet  problem  which 
specifies  a  flush-type  inlet,  it  might  be  anticipated  that  non- 
uniform  inlet  velocity  profiles  could  be  generated  with 
corresponding  diffuser  performance  problems.  This  is  only 
one  example  of  the  type  of  feature  which  a  given  geometry 
might  impose  on  the  flowfield,  a  result  which  becomes  clear 
and  casual  in  this  inverse  formulation. 
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Inverse  Design  of  Optimal 
Diffusers  With  Experimental 
Corroboration 

TI:e  equations  of  motion  for  a  turbulent  boundary  layer  are  formulated  in  an  in¬ 
verse  manner  and  solved  to  provide  the  geometrical  configuration  of  optimal  two- 
dimensional  diffusers  with  a  distinct  inviscid  core.  For  each  Reynolds  number/ dif¬ 
fusion  length,  optimality  is  inferred  when  the  lowest  exit  freestream  velocity  is  com¬ 
puted  f>  om  the  members  of  an  assumed  one-parameter  family  of  skin  friction  decay 
functions.  In  addition  to  the  unique  inverse  formulation  and  establishment  of  the 
appropriate  skin  friction  distribution,  an  integral  part  of  the  analysis  is  the  develop¬ 
ment  of  an  invertable  adverse  pressure  gradient  skin  friction  law,  which  is  used  to 
solve  for  the  channel-averaged  streamwise  pressure  ( inviscid  core  velocity) 
distribution.  The  diffuser  configurations  are  presented  graphically  for  various 
Reynolds  numbers/ lengths/ blockages  and  exhibit  a  distinctive  concave  shape.  The 
results  are  compared  to  three  sets  of  experimental  data  which  support  the  validity  of 
the  viscous  flow  analysis,  the  formulation  of  the  governing  equations,  and  the  skin 
friction  law  which  is  proposed. 


Introduction 


Theoretical  determination  of  the  wall  contour,  which  is 
required  to  establish  a  specified  flowfield  diffusion  and/or 
turning,  is  an  elegant  application  of  the  inverse  potential  flow 
theory  in  (1,  21.  In  this  approach  the  surface  pressure 
distribution  in  the  transformed  mathematical  plane  miu  be 
input  to  the  analysis.  The  so’ution  and  subsequent  retrans¬ 
formation  to  the  physical  plane  provides  the  desired  duct 
configuration.  In  order  to  specify  a  meaningful  surface 
pressure  distribution,  however,  the  boundary  layer  behavior 
must  be  considered. 

In  recent  years,  turbulence  modeling  methodology  (see  (3)) 
has  been  elevated  to  the  point  where  equations  which 
represent  the  entire  Reynolds  stress  tensor  may  be  applied  to 
compute  the  flow  through  any  specified  channel.  Ad¬ 
ditionally,  available  experimental  data  [4]  on  diffuser  per¬ 
formance  is  coming  »nder  control  with  the  proper  treatment 
of  boundary  layer  e.i'ects.  With  such  correlated  data,  the 
general  configuration  'f  optimal  diffusers  may  be  inferred. 
Ideally,  the  most  accurate  analytical  models  should  be  utilized 
to  compute  the  detailed  diffuser  geometry.  Complexities  in 
such  an  analytical  formulation  and  its  computational 
magnitude  have,  to  this  point,  prohibited  such  an  approach. 
In  fact,  few  analyses  are  available  which  provide  a  wholly 
rational  diffuser  design  and  those  which  seek  out  the  optimum 
configurations  are  rare. 

Generally,  an  unseparated  flow  condition  with  maximum 
static  pressure  recovery  in  the  diffuser  is  sought.  The 
displacement  thickness  distribution  is  computed  to  correct  the 
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indicated  pote.i  ial  flow  wall  contour  for  boundary  layer 
blockage  effects.  In  the  direct  problem  this  procedure  is 
iterative  (not  necessarily  convergent)  and  provides  no 
straightforward  mechanism  for  establishing  an  optimal 
diffuser  contour. 

To  fill  this  void,  a  new  diffuser  design  technique  has  been 
developed  •*  hi  eh  utilizes  an  original  algebraic  skin  friction 
model  based  on  data  correlations  for  flows  in  adverre  pressure 
gradients.  The  conventional  boundary  layer  equations  are 
combined  with  this  skin  friction  model  in  ah  inverse  manner 
and  optimality  is  established  by  maximizing  the  diffuser 
pressure  recovery  with  respect  to  a  single  parameter  which 
characterizes  the  skin  friction  distribution  along  the  diffuser. 
The  governing  equations  then  provide  simultaneously  the 
boundary  layer  thickness  and  freestream  velocity  distribution, 
from  which  the  optimal  diffuser  configuration  may  be  in¬ 
ferred  from  either  one-dimensional  or  inverse  potential  flow. 

In  conjunction  with  available  experimental  data,  the 
analysis  establishes  a  rational  diffuser  design  methodology, 
affirms  the  validity  of  a  new  skin  friction  corrcla'ion,  defines 
the  characteristic  optimal  skin  friction  decay,  and  generates 
conceptually  new  diffuser  wall  contours  for  experimental 
evaluation. 

Analysts 

Since  an  unseparated  diffuser  is  sought,  (he  skin  friction 
should  vary  from  its  constant  pressure  value  near  the  diffuser 
entrance  to  a  small  value  at  the  exit.  The  optimal  path  bet¬ 
ween  these  two  extremes  is  that  which  provides  the  lowest 
inviscid  core  velocity  at  the  diffuser  exit  (i.*„  maximum  static 
pressure  recovery).  In  order  to  restrict  the  mathematical  scope 
of  the  problem,  the  optimization  process  is  simplified  by 
specifying  tite  skin  friction  decay  function  to  be  a  one 
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parameter  family  which  approaches  zero  with  a  zero 
derivative  at  the  diffuser  exit.  This  exit  condition  enhances  the 
possibility  of  establishing  an  equilibrium  incipiently  separated 
flow  as  in  [51,  but  excludes  cases  of  increased  recovery  with 
slight  stall  in  (6]  which  have  their  origin  in  explicit  transient 
flow  phen  OTena. 

Such  a  one  parameter  family  with  the  requisite  boundary 
conditions  is  given  by  the  equation. 

r,=(i-i)'  in 

By  assuming  a  distribution  of  values  for  A'  ( 2r  1 )  lor  each 
flowfield  studied,  the  equations  of  motion  may  be  solved  to 
determine  simultaneously  the  corresponding  distributions  of 
&’  and  U  along  the  channel.  The  A'  value  which  provides  the 
minimum  freestream  velocity  at  the  channel  exit  is  taken  to  be 
that  which  characterizes  the  optimal  skin  friction  decay.  It 
will  become  apparent  that  the  optimal  A'  is  not  infinitely  large, 
and  this  fact  can  be  attributed  primarily  to  the  shape  factor 
growth  as  separation  is  approached. 

Clearly,  the  assumed  skin  friction  distribution  must  be 
input  to  the  equations  of  motion  which  are  utilized  here  in  the 
following  integral  form: 

0(a)  -  0,  =  ( *  <  C,Cj_ /2  -  S( H -v 2 )S,,  )dx  (2) 

Note  that  the  pressure  gradient  term  is  included  in  the  in¬ 
tegrand  on  the  right-hand  side  of  equation  (2)  due  to  the 
inverse  nature  of  the  computational  procedure,  as  discussed  in 
the  next  section.  Additionally,  expressions  for  H  lin  lerms  of 
H.j.H,-,.  and  5)  and  Slin  terms  of  C,  i  must  be  developed.  The 
nonlinear  interpolation  formula  developed  in  [7]  for  H  was 
utilized,  with  values  for  H ,,  hc ,  and  S-,  as  functions  of 
Reynolds  number  taken  from  [7]  through  (13),  along  with  the 
C,  correlation.  By  correlating  the  data  from  (7),  the 
following  adverse  pressure  gradient  skin  friction  law,  which 
forms  the  basis  for  all  the  analytical  results  that  follow,  is 
obtained: 


e,= 

1-5 

l  —  ■  '"■! 

uoi 

i-6 

(3) 

i  *  06 

i 

)M  -- 

j-  2g(Re„) 

(d) 

g*(Ret,r 

3 »  0 . 70:  -  0. 1  os  (log,,  Re, 4 

10  Re„):  (5) 


The  (enforced)  asymptote  for  <t>  is  10.22  at  Re„  =  5.76  x  10s, 
which  is  the  high  Reynolds  number  (limiting  flow)  value 
appropriate  to  the  developments  in  [7] . 

Future  work  utilizing  algebraic  skin  friction  correlations 
which  contain  explicit  effects  due  to  turbulent  flow  history 
and  wall  curvature  should  increase  the  universality  of  this 
analytical  approach. 

Method  of  Solution 

In  order  to  integrate  equation  (2)  with  conventional 
quadrature  algorithms,  the  characteristic  flowfield 
parameters  ReH  and  Ll 9,  must  be  specified,  allowing  the 
constant  of  integration  (0, )  and  to  be  determined.  Ad¬ 
ditionally,  equation  (I)  provides  C,  once  A’  is  specified.  The 
heart  of  the  computation,  however,  lies  in  utilizing  equation 
(3)  to  supply  S  in  equation  (2),  in  terms  of  the  known  Ct.  With 
this  mating  of  equations  (1)  through  (3)  and  the  functional 
formulation  for  H,  equation  (2)  may  be  integrated  utilizing 
expressions  for  C,„,  H,„  H„  and  Ser  in  terms  of  the  local 
Reynolds  number  (Re„).  In  order  to  determine  the  local 
Reynolds  number,  the  local  freestream  velocity  must  be  , 
known  and  is  obtained  from  the  concurrent  integration  of 
equation  (3),  utilizing  the  aforementioned  functional  inputs. 
The  treestream  velocity  is  also  required  to  establish  the  dif¬ 
fuser  geometry  based  on  an  inviscid  flow  analysis. 

From  equations  (1)  and  (3),  and  the  definition  of  &  it  is 
apparent  that  rapid  growth  of  6  is  detrimental  to  obtaining 
large  freestream  velocity  decelerations.  The  optimal  0  growth 
is  a  trade-off  between  the  rate  at  which  5  approaches  one  (i.e., 
Ct  —  0)  in  equations  (2)  am!  (3).  and  the  rate  at  which  H  in¬ 
creases  to  in  equation  (2).  Therefore,  an  optimal,  finite 
value  of  A’  is  to  be  expected. 

In  summary,  after  specifying  Re,,,  and  L/B  .  the  solution 
of  the  governing  equations  provides  the  distributions  of  0,  6, 
and  H  for  a  given  A.  Upon  solving  for  0,  as  a  function  of  A‘ 
at  a  given  Re(,  -  L/B,  combination,  the  optimal  diffusion 
(for  this  particular  one  parameter  family  of  skin  friction 
decay  functions)  may  be  identified.  In  order  to  display  the 
geometric  configuration  of  these  optimal  diffusers,  C'  is 
assumed  to  be  uniform  across  the  inviscid  core  and  the  one- 
dimensional  continuity  equation  is  applied  to  a  two- 
dimensional  diffuser  (divergence  in  one  plane)  resulting  in: 

WimK/L)  4-UB  i*  (6) 

where  the  physical  imerpretaiion  of  displacement  thtekness 
has  been  utilized.'  In  addition,  the  effective  inviscid  core 
height  at  the  inlet  has  been  utilized  as  a  characteristic  length 
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Nomenclature 


AH  » 
B  - 
B  - 

C.  » 

6  * 
O',  » 
c .  *• 

i  • 
H  « 
h  « 


area  ratio 
(A*  ft) 

vkm  friction  coefficient 

C,»C. 

CMC,. 

diffuser  static  p?c>vwte  recovery 
coefficient 

funct  or  of  He.,  in  equation  t.5> 
shape  idcior.  i*  O 
effective  inviscid  core  half- 
height 


L  *  Mreamwise  wetted  length 

V  «  exponent  inequation  (I > 
Re  V  Reynold*  number 

5  «  *>■  d  ir.  C  t  d\ 
i  *  fcS„ 

V  *>  frertream  vciocH) 

C  ct  (  '  { •' 

u  a  channel  half  heiRlu 
a  >t rcjtJiiw  i m;  coordinate 

x  «  x  L 

h*  v  dupiacerncm  thickness 

i'*  *  »*  i; 
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0  «t  momentum  thickness 

0  »  9  0, 

<s  »  function  of  gin  equal  ton  (4) 
SuHH'rtpiS 

er  *  condition-  at  separation 
/  «  conditions  at  diffuser  exit 
i  *  conditions  s'  diliuvct  mlci 
os  constant  pressure 
0  «  referenced  to  0  ■ 

L,  *  referenced  to  upstream  wetted 
length 
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(S'  »  (4V*>  ).  in  !icu  of  (he  total  channel  height,  to  obtain  a 
concise  functional  form  with  linear  dependence  upon  a  single 
boundary  layer  parameier.  The  conventional  boundary  layer 
blockage  (£  «($*•'»),)  is  related  10  B'  by  Bo  BVtl  ♦  B’l 
By  specifying  ;he  diffuser  blockage,  or  in  viscid  core  as¬ 
pect  ratio  {Lfhj  which  are  related  by  L/h.  * 
(I/O,  MB’///, (Re.,  )  |.  equation  (6)  was  evaluated  with  the 
solutions  of  the  equations  of  motion  to  provide  the  channel 
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Ftg.  4  Diffuser  will  and  Invlscld  cor*  art*  distribution  for  R#8, 
*105-10‘ ati/S|  »104  and#'  *0.01,  U2W,  -3770,39.41,40.71. 
41.71 


Ftg.  5  OltfuMK  wall  and  Inviacid  cor*  area  distribution  tor  Roe, 
-10*  -10*  at  k/t>(  >10*  and  r  >0.0011/21*,  >  79.01.79.41 9107. 
94.09 


geometry.  The  hi  viscid  core  aspect  ratio  is  related  to  the 
conventional  diffuser  aspect  ratio  (£/2M',)  by  L2H‘, 
■  (1/A,  I  (1 4*fl‘) ' ‘<*2  (Z./A.H1  -  BV2.  Note  that  other 

commonly  used  variables,  such  as  duct  Reynolds  number,  are 
simple  functions  of  the  parameters  utilized  here  to  charac¬ 
terize  the  Ho*  field. 

When  the  diffuser  geometries  are  presented  graphically,  the 
B‘  ■  0  curve  ts  plotted  to  represent  the  effective  invisetd  core 
area  distribution.  A  second  curve,  which  is  the  t assumed) 
invisetd  core  added  to  a  linearly  scaled  and  normalized 
displacement  thickness  distribution,  is  plotied  to  represent  the 
(two-dimensional)  diffuser  wall,  therefore,  with  this 
boundary  layer  blockage  representation,  the  displacement 
thickness  growth  is  readily  observed  and  rescaled  for  different 
flow  Held  conditions. 

Clearly,  the  computed  pressure  gradient  could  be  estab¬ 
lished  wiihout  a  distinct  invisetd  structure  and  the  resultant 
diffuser  area  gradicnts/divergcnce  angles  would  be  smaller. 
In  such  computations,  it  is  more  relevant  to  associate  the 
definition  of  $  directly  with  the  pressure  gradient  and  utilize 
the  channel-averaged  momentum  and  continuity  equation*  to 
provide  the  cross-sectional  area  distribution  for  any  geometric 
class  of  diffusers. 

fn  either  case,  the  one-dimensional  flow  assumption  carries 
the  inherent  approximation  of  neglecting  flow  non- 
uniformities  across  ihe  diffuser  due  to  wall  curvature.  This 
effect  shall  be  assessed  in  turure  research  utilizing  inverse 
potential  flow  analysis  (l,  2)  to  generate  the  wall  contour. 
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Analytical  Results 

For  Reynolds  numbers  from  Re0  =  10'  -10"  and  dif¬ 
fusion  lengths  from  i/0,  =  10'-10\  the  governing 
equations  were  solved  numerically  utilizing  a  distribution  of 
values  for  S'  sufficient  to  determine  the  optima)  skin  friction 
decay  in  each  situation.  The  velocity  ratio  across  the  diffuser 
(for  the  optimal  S)  has  been  plotted  versus  Rew  with  i/0, 
as  a  parameter  in  Fig.  1.  For  a  given  diffuser  length  and  inlet 
Reynolds  number,  reducing  the  diffuser  entrance  momentum 
thickness  is  a  process  which  moves  diagonally  across  the 
figure  toward  the  lower  left  hand  corner. 

The  diffuser  geometries  corresponding  to  the  decelerations 
indicated  in  Fig.  1  are  presented  in  Figs.  2  through  5  for  i/0, 
=  10:  - 10'  and  Re(lj  =  10'  -  I0k  which  encompasses  the 
configurations  of  primary  interest.  Note  that  each  relative 
length  configuration  is  associated  with  a  specified  blockage 
factor  (0'),  which  was  selected  for  clarity  of  data  presen¬ 
tation.  The  corresponding  aspect  ratios  (as  indicated  on  the 
figure  caption!  are  not  necessarily  always  physically  realistic 
and  the  data  may  be  re-scaled  to  any  desired  blockage. 

Of  particular  interest  is  the  fact  that  the  diffusers  which 
have  been  computed  assume  a  characteristic  concave  shape 
(when  viewing  the  flow  side,  excepting  ihe  immediate  en¬ 
trance)  so  as  to  retard  the  thin  upstream  boundary  layers  more 
strongly  than  the  downstream  flow.  This  is  an  extension  of  the 
logical  deceleration  mechanism  which  has  been  set  forth  in 
proposing  and  rationalizing  the  performance  of  straight-wall 
diffusers.  The  unique  diffuser  configuration  at  the  entrance 
plane  is  a  result  of  the  lag  between  constant  pressure  initial 
conditions  and  the  establishment  of  a  boundary  layer  growth 
in  response  to  strong  adverse  pressure  gradients,  Ex¬ 
perimentally  observed  conditions  (5]  at  the  immediate  en¬ 
trance  to  the  diffuser  are  also  associated  with  the  transition  of 
the  pressure  field  from  constant  to  increasing  pressure. 

The  sensiriviry  of  the  results  to  Ain  the  neighborhood  of  its 
optimum  value  may  be  seen  in  the  data  of  Ftg,  6.  Clearly,  the 
boundary  layer  growth  is  sensitive  to  the  pressure  gradient 
but,  utilizing  the  inverse  formulation,  the  optimal  pressure 
gradient  is  not  sensitive  to  .V  in  the  netfhhorhond  of  the 
optimum  S.  Therefore,  the  results  presented  here  should  not 
be  unduly  restrictive  due  to  the  choice  of  equation  ( I ). 

Experimental  Comparisons  and  Conclusions 

The  analytical  revulis  presented  here  provide  a  graphic 
indication  of  the  degree  to  which  boundary  layer  control 
allows  even  targe  area  ratio  diffusers  to  perform  efficiently. 
Consider  tbe  dramatic  results  from  |15|  in  which  the  static 
pressure  recovery  efficiency  for  a4:l  area  ratio  diffuser  sarin 
from  ©.?0  at  0  «  0.06  to  approximately  0.0?  with  uniform 
entrance  flows  (inferred  by  extrapolation).  Maintaining  a 
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turbulent  entrance  boundary  layer  and  an  incompressible 
flow,  the  diffusion  length  for  these  blockages  covers  the  range 
1/0,  =  3  x  102  — 105  at  a  Reynolds  number  of  10’.  The 
results  of  the  comnutations  carried  out  here  indicate  that  the 
efficiency  of  an  optimal  7.15:1  aspect  ratio  <1/2  W,)  diffuser 
(as  used  in  (15J)  increases  from  0.75  to  0.99  over  this  variation 
in  1/0,  as  the  entrance  momentum  thickness  essentially 
approaches  zero.  This  corroboration  supports  the  position 
that  very  large  optimal  decelerations  predicted  theoretically 
are  experimentally  achievable  for  sufficiently  thin  diffuser 
entrance  boundary  layers. 

The  extensively  quoted  data  of  |16J  provide  a  further  test  of 
the  analytical  results  obtained  here.  For  a  blockage  factor  (0) 
of  0.03  and  an  aspect  ratio  (L/2W, )  of  18.0,  the  experimental 
results  in  Fig.  10  (16)  indicate  an  optimum  C,  of  0.73  at  an 
area  ratio  of  3.0  which  leads  to  an  efficiency  of  0.82.  The 
value  of  Retli  appropriate  to  the  experiment  in  (13)  is  ap¬ 
proximately  1.5  x  10J  with  a  corresponding  L/Q,  -  1.8  x 
10’.  If  the  data  obtained  here  at  Re0  ■  10'  andZ./0,  »  IQ1 
is  utilized,  the  optimum  diffuser  is  indicated  at  an  area  ratio 
of  3.25  with  an  efficiency  of  0.82.  The  theoretical  results  are, 
therefore,  reasonable  predictions  of  the  actual  diffuser 
performance  and  reinforce  the  experimental  indications  from 
(16)  of  a  potential  benefit  from  contouring  the  skin  friction 
decay  and  resultant  diffuser  cross-sectional  area  distribution. 

A  targe  portion  of  the  experimental  studies  on  diffusers 
reported  in  the  open  literature  lie*  in  the  regime  L/Q,  »  I0! 
-  10*  and  Reft(  ■  10’  -  10*.  The  optimal  two-dimensional 
diffusers  which  have  been  plotted  look  ’•familiar*’  for  the 
relatively  larger  Q,  values  In  this  range.  In  fact,  as  has  been 
shown,  the  analytically  computed  diffuser  performance 
closely  follow*  available  experimental  results.  It  is  apparent 
that,  as  the  design  point  for  the  diffuser  penetrates  the  low 
Rcw,  •  high  L/Q,  corner  of  the  analysis,  the  configurations 
depart  from  those  commonly  accepted  for  optimal  diffusers. 
Therefore,  the  excellent  NBS  data  from  (5)  has  been  used  to 
validate  the  analysis  for  large  diffusion  lengths  and 
freestream  velocity  decelerations,  inasmuch  as  it  foiges 
further  into  the  clean  entrance  domain  with  essentially  two- 
dimensional  (low.  That  experimental  work  was  also  directed 
at  establishing  an  incipiemly  separating,  yet  stable,  boundary- 
layer  flow  which  is  relevant  to  the  class  of  optimal  diffusers 
presented  here. 

From  the  wind  tunnel  turbulence  damping  screen  to  the 
streamwise  position  where  (m  Fig.  3  of  (5|)  the  adverse 
pressure  gradient  takes  hold  is  approximately  1.45?  m.  Based 
on  the  conditions  quoted  in  (5).  the  resultant  value  of  0, 
relevant  to  the  analysis  employed  here  is  2.881  mm,  which  was 
obtained  from  the  flat  plate  turbi  lent  boundarv  layer  growth 
’’law”:  0,/i,  «  0,036/Re*  ej.  This  agrees  with  the  ex- 
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Fig.  5  Compart  son  of  analytical  and  axpadmanfal  displacement 
thickness  distributions  at  l/#,  •  2.787  *  10*  and  Refy  a  3.955  x  103 


pcrimental  results  in  Fig.  8,  (5]  and  establishes  a 
corresponding  diffuser  entrance  Reynolds  number  of  Re#  = 
3.955  x  103. 

With  the  exception  of  the  transition  from  constant  pressure 
to  an  adverse  pressure  gradient,  the  skin  friction  coefficient 
profile  from  Fig.  5  {5]  is  similar  to  the  decay  function  which 
was  utilized  here.  It  is  to  be  remembered  that  the  dependence 
of  Of  (and  the  resultant  streamwise  profile)  on  N  is  relatively 
weak  in  the  neighborhood  of  the  optimum  <V.  Therefore,  the 
experimental  reality  is  that  the  skin  friction  decay  in  (5(  which 
establishes  incipient  separation  may  be  optimal. 

The  results  in  Fig.  5  (51  indicate  that  the  point  of  zero  skin 
friction  would  occur  at  approximately  8.13  m.  This  implies  an 
4/0,  *  2.757  x  t0J  after  accounting  for  the  0.187-m 
transition  section  between  the  uniform  velocity  and  power  law 
velocity  decay  regions.  Utilizing  an  Rc^  -  3.955  x  I03  and 
4/e,  ■  2.757  x  10\  the  data  from  Fig.  1  may  be  entered 
with  a  non-linear  interpolation  scheme  which  provides  a  value 
for  0/  m  0.?9  *  0.01.  The  results  of  Fig.  3  (5|  may  be  ex¬ 
trapolated  to  8.13  m,  indicating  a  0,  «  0.38  *  0.02,  In  fact, 
the  entire  freestreem  velocity  and  displacement  thickness 
distributions  predicted  by  analysis,  shown  in  Figs.  7  and  8, 
compare  favorably  with  that  which  (5J  presents  as  charac¬ 
teristic  of  incipient  separation. 

In  summary,  the  analytical  results  obtained  here  are 
substantiated  by  experimental  data  for  flows  in  adverse 
pressure  gradients  which  approach  separation  and,  therefore. 


in,  miuiij  oi  ct|uduuiib  v  i  j  ana  pTaim  me  overall  mcoiciicai 
viscous  flow  approach  is  confirmed.  The  geometrical  con¬ 
figuration  of  the  diffusers  obtained  from  quasi-one- 
dimensional  flow  analysis  points  to  the  need  for  a  critical 
assessment  of  diffuser  contouring  (as  in  [16])  and  its  in¬ 
teraction  with  the  incoming  boundary  layer.  In  addition,  the 
examination  of  large  diffusion  in  the  presence  of  very  clean 
entrance  conditions  (B  -0(10  "*))  through  use  of  efficiently 
designed  bundary  layer  cut-off  slots  [2]  is  indicated. 
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Nondskw 
■speed  of  sound 
■cross-sectional  area 
■constant  in  Eq.  (7)  for  4 

■  particle  drag  coefficient 
■skin  friction  cocfficicat 
■particle  scale 
■hydraulic  diameter 
■Euler  number,  p/pa*1 

«  functions  defined  in  Eq.  (7) 
■generalized  Integrand 
■constant  in  Eq.  (I),  2  Ci(UD%) 
■constant  in  Eq.  (I),  3  CD(L/Dt)/4 

■  nozzle  length 
■Mach  number 


■pressure 

■  ratio  of  liquid  to  gas  mass  fiow  rate 

■  nozzle  surface  area  per  unit  length 
■veiodty 


■relative  velocity,  u(  -a, 
■bulk  velocity 
•  volume  ratio,  (p,/pJ/r 

■  Weber  number,  pau>t/« 

■  stream  wisepoeition 
-A#* -I 

■parameters  defined  by  Eq.  (7) 

■  (p(a;/p ,uf)* 


■density 
■surface  tension 
■momentum  thickness 


■efficiency.  (V<Uw)1 


Marrjm 


oil  ■conditions  for  partide  breakup 

e  «  nozzle  out 

?■  gaseous 
■liquid 


0  ■  initial  or  reference  conditions 

t  ■stagnation  conditions 

Superscripts 

(  )*  ■critical  or  throat  conditions 
(-)  ■normalized  by  throat  conditions  or  L 

(')  ■differentiation  with  respect  to  x 

la  trod  taction 

COMPLICATIONS  inuoduced  by  multiphase  flow  to  the 
design  of  efficient  converging-diverging  nozzle*  are 
significant. 'A  Optimal  length  scaling  (aspect  ratio),  coo* 
touring,  and  area  ratio  arc  sensitive  to  the  detailed  bubble 
growth/ partick  breakup  and  dynamics.  As  discussed  in  Ref. 
3,  these  phenomena  are  often  complea  and  not  easily 
modeled.  This  restricts  the  utility  of  large-scale  computational 
algorithms  whose  design  predictions  art  difficult  to  distill  into 
a  compact  form.  Therefore,  an  approximate  solution  which 
perturb*  the  iseatropic,  homogeneous,  equilibrium  (IHE) 
flow  model  was  developed  to  synthesize  the  optimal  nozzle 
geometry  for  anisentropic.  tummy-pus,  nooequillbeium 

situations 

Area  Ratio 

In  the  limit  of  large  volume  ratio  (it,),  and  assuming  that 
the  density  ratio  is  sufficiently  small  for  the  mass  ratio  to  be 
>  I,  the  (HE  model  of  two-phase  now  (ice  Ref.  4)  may  be 
expanded  to  provide  the  following  approximate  relationships: 


«tt,**(p(/p) 

Pi 

(1) 

40d 

0p‘tl 

~Z —  -v,/2 

Pt 

(2) 

Applying  Eq.  (1)  at  the  nozzle  throat  and  substituting  for  a 
from  Eq.  (2)  results  in 

(1) 
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Utilizing  continuity  of  mass  flow,  the  local  area  ratio  becomes 

A}=e*/U+P)  (5) 

where  fimU(p,/p)}  -\  and,  therefore,  from  Eq.  (4) 
Q=M3  - 1. 

The  validity  of  these  formulas  may  be  evaluated  by  com¬ 
parison  with  the  results  in  Refs.  1  and  2  fot  air/water  nozzles 
operating  at  r *  10  and  p,lpt* 35.  Equation  (5)  predicts  full 
expansion  at  ^,=7.96,  whereas  the  (suboptimlfl)  pressure 
contour  design  analysis  in  Ref.  2  results  in  A,* 7.022. 
Evaluation  mode  results  were  computed  in  Ref.  1  for  identical 
conditions  in  an  annular  nozzle  with  an  approximately  linear 
area  distribution  in  the  supersonic  region.  These  com¬ 
putations  indicate  that  sensible  expansion  ceases  at  A  -6.93 
for  both  7  and  21  deg  divergence  angle  nozzles.  Equation  (3) 
satisfactorily  predicts  the  throat  conditions  for  both  sets  of 
data. 

This  simple  IHE  limit  model  provides  the  framework  within 
which  anisentropic,  nonhomogeneous,  nonequilibrium  flow 
phenomena  may  be  analyzed.  Based  on  the  substantiation  of 
Eqs.  (3)  and  (5),  analytical  representations  of  the  two-phase- 
flowfield  properties  allow  for  the  evaluation  of  two  “loss” 
mechanisms  of  interest:  nozzle  wall  friction  and  interphase 
drag. 

It  will  be  shown  that  optimal  contouring  based  on  wall 
shear  minimization  is  derived  as  a  function  of  nozzle  length 
scale.  In  addition,  the  essence  of  the  interphase  dynamics 
analysis  is  to  optimize  the  nozzle  length  scale  in  a  tradeoff 
between  wall  friction  and  phase  slip.  Inasmuch  as  the  optimal 
nozzle  contour  (combined  with  the  IHE  pressure-area 
relationship)  is  used  to  provide  the  bulk  flow  pressure 
distribution  for  the  particle  dynamics  analysis,  there  is  a 
coupling  between  the  following  contouring  and  length  scale 
estimations.  Therefore,  although  the  contouring  (via  wall 
shear  minimization  considerations)  and  length  optimization 
(via  exit  bulk  velocity  maximization— Including  wall  shear, 
interphase  drag,  and  pressure  effects)  of  the  nozzle  are  carried 
out  separately,  they  are  not  independent  and  the  final  optimal 
configuration  reflects  the  integrated  effect  of  all  phenomena. 


Fig.  I  /mi  My*  A. 


Fig. 2  Hp,/p)mik\JinA. 


Contouring 

In  order  to  establish  a  nozzle  contour  which  “encourages” 
the  IHE  flowfield,  an  optimization  procedure  was  developed 
in  which  the  wall  friction,  Integrated  over  the  nozzle  surface 
area,  Is  minimized.  This  procedure  may  be  compared  to  the 
approach  in  Ref.  2  wherein  a  specific  family  of  nozzle 
pressure  profiles  is  select td  a  priori.  The  particular  profile 
which  provides  the  maximum  exit  bulk  velocity  is  assumed  to 
be  optimal.  These  approaches  ere  also  to  be  contrasted  with 
classical  wave  cancellation  contouring  of  gas  flow  nozzles. 

Applying  traditional  variational  principles  to  the  integral 
Jew'CydS,  the  following  equation  for  the  eroei -sectional  area 
distribution  of  nozzles  (wheee  shape  remains  similar 
throughout)  may  be  derived 
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In  order  to  derive  Eq.  (6),  the  Eukr-Lagraagc  equation: 
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which  maternal  the  minimum  value  for  a  definite  integral 
is  tppUad  to  bu*Cjd&  In  this  case,  /■>£" 
(l+h)  where/laciadm  the  dynamic  pressure  and  local  body 
radius  Cffsct,  g  accounts  fot  the  change  la  local  wetted  area 
due  to  the  of  the  aozzle  surface,  and  A  accounts  for 
changes  ta  C.  dot  to  pressure  gradients.  Note  that  Eqs.  (1*9) 
allow /and  /  to  be  expremad  as  functions  of  A  and  Eqs.  (7) 
convert /ton  function  of  A  and  A  with  sand  das  panunsun 
In  the  optimization.  These  characterise  the  nozzle  aspect  ratio 
and  sansttivity  of  wall  friction  to  pressure  gradient,  respec¬ 
tively. 

Figure  I  pretests  both  subsonic  and  supersonic  branches  of 
/and  U  as  a  function  of  arm  ratio,  with  the  maximum  value 
of/ocotrriag  at  A •  1.43  and  Af* 1.74.  The  corresponding 
results  foe  /  and  p,tp  are  presented  in  Fig.  2.  Given  an  op¬ 
timal  pressure  distribution,1  the  nozzle  contour  Is.  therefore, 
known.  However,  ia  this  research,  the  optimal  area 
distribution  (for  a  given  length  scale)  shall  be  iafemd  from 
the  titbit  ioa  of  Eq.  (4)  and  the  resultant  implied  pressure 
gradient  used  to  optinsiit  the  length  scale  based  oa  subsequent 
phase dip  comidcratioas. 


and  Phase  Dynamics 

The  IHE  Emit  model  and  octbca)  nozzle  arm 

*m(4UDt)-,t  4<ac £/(i/P*)(  /«(d//d4  )/(#**)  <»)  distribution  allow  the  panick/gas  dynamics  to  be  approached 
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la  a  particularly  con  veal  eat  way,  to  order  to  complete  the 
optimization  peocedur*.  Under  ihemtamptiou  of  high  e„  the 
IHB  pmturMiM  relationship  tad  optimal  mn  distribution 
ut  input  10  the  analysis  to  provide  pressure  ud  am  h  ■ 
taction  of  £  for  a  given  amklsagta  scale  (aspect  ratio),  la 
this  com,  the  phase  ittp  tad  Uqnld  veiodtias  an  coopwad 
from14 


hjd£/df 

Mf^+x,*,) 


•dtM  joint  {fetation  with  da 

optisnumann distribution  (inquired.  - 


At  this  point,  validation  -af  the  approach  reties  prisasrily  oo 
the  data  lo  Ref.  2,  wheraia  two  different  air/wnttr  amoks 


w«e  analyzed  and  tested  at  (p,/pt)  >33  with  r>  10.  These 
nozzles  had  identical  subsonic  contraction  cones  and  two 
different  length  supersonic  sections  (7  vs  21  deg  divergence 
angles),  with  the  same  overall  area  ratio.  Therefore,  the  area 
contour  and  length  optimization  developed  here  will  also 
focus  on  the  supersonic  regime.  This  avoids  complex  particle 
breakup  dynamic*  modeling,  inasmuch  as  liquid  particulates 
typically  reach  a  limiting  *l*e  at  the  throat. ,J 

Given  a  reliable  particle  breakup  model  (e.g„  *>> 
the  entite  nozde  may  be  computed  with  the  same  contour  and 
particle  dynamics  analysis  (utilizing  dm  local  particle  scale), 
Inasmuch  as  the  IHB  pressure-area  relationship  is  applicable 
to  the  entire  nozzle.  For  cases  of  fUed-si»  particles,  the  entire 
«««h  my  hg  csunutad  mihmws  these  Muatiottt  without 
modification. 

Typical  integration  of  Eq.  (g)  from  the  throat  to  the  exit 
plane  is  presented  in  Fig.  3  and  indkatss  that  a  simple  (nearly) 
linear  area  distribution  is  often  optimal.  (For  axhymmetric 
nozzles,  tide  results  in  a  configuration  approaching  a  ball 
shape.)  Tba  similarity  between  these  shapat  and  thorn  in  Ref. 
2  is  remarkable.  This,  combined  with  the  substantiation  of 
Eqs.  (3)  and  (f)  indicates  that  the  information  provided  in 
Ref.  2  on  the  effect  of  nozzle  length  should  also  be  relevant  to 
the  tafeth  scaling  optimisation  procedure  developed  here. 

Note  that  the  rotation  to  6q.  (»)  is  •  spirt-point  boundary 
value  problem.  A*  •  1  it  required  by  ostutition  of  A  and  the 
origin  of  integration  at  the  throet.  Either  the  subsonic  or 
supettoric  «|Im«  n,«y  ge  evaluated,  depending  on 
wNch  branch  of /and /era  chosen  tithe  integration  proceeds 
awey  from  the  thraet.  Foe  a  given  i  end  I,  in  order  to  obtain 
the  desired  A,,  the  appropriate  4,  is  determined  iteratively. 
In  lieu  cf  Util  iteration,  the  computations  may  be  carried  out 
whhvig  as  a  parameter.  A  family  of  nozzle  ihapas/ane  ratios 

U  wuvwyjHM* 

For  e  (p,/p()>3 3  end  r>IO,  the  IHB  pressure«ren 
raietionshipe  era  combined  with  the  optimal  contour  to  solve 
Eqs.  (f)  and  gtaerau  profiles  of  liquid  end  slip  velocities 
along  the  node  for  various  values  of  Ar,  and  k,  typical 
results  are  shown  In  Fig.  4,  whkh  Uhutnte  a  dependence  on 
ftotde  length  slmilec  to  that  prntniid  In  Rtf.  7  The  fugitive 
vahtm  of  a,  shown  for  Ar,  «•  IG*  and  reflect  e 

somewhat  nonphysical  disparity  between  C>  end  C0,  with  the 
well  sheer  tueti  dominating. 

It  may  be  shown  that  the  hulk  (taase  flow.wdghtad  mean) 
velocity  fat  rdated  to  the  slipaad  liquid  velocities  by 

<i  •«,-<*,/(/+')  (10) 
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The  measured  bulk  velocities  in  Ref.  2  are  170  and  183  m/s 
for  nozzles  I  (7  deg)  and  I!  (21  deg)  and  resulted  in  efficiencies 
(a)  of  0.48  and  0.56,  respectively.  Utilizing  Eqs.  (9)  and  (10), 
the  computed  nozzle  efficiency  is  plotted  in  Fig.  5  as  a  func¬ 
tion  of  k,  for  various  k}.  Compared  with  the  data  from  Ref. 
2,  it  is  seen  that  the  trend  toward  higher  efficiency  with 
shorter  nozzles  is  predicted  analytically.  For  these  flow 
conditions,  this  dependence  is  observed  only  for  large  k, 
nozzles,  due  here  to  the  small  throat  hydraulic  diameter  of 
annular  nozzles. 

Conclusion 

An  approximate,  self-consistent  procedure  has  been 
developed  to  determine  the  optimal  area  ratio,  contour,  and 
aspect  ratio  of  two-phase-flow  nozzles  with  large  volume 
ratios.  Results  of  the  analysis  have  been  validated  by  com¬ 
parison  with  experimental  data  obtained  from  nozzles  which 
closely  approximate  the  computed  optimal  configuration.  At 
this  stage  of  the  analysis,  fixed-size  particles  were  considered 
in  the  assessment  of  optimal  nozzle  aspect  ratio.  However, 


given  a  reliable  particle  size  model  for  any  given  situation,  the 
identical  procedure  may  be  employed  utilizing  a  variable 
particle  scale. 

Acknowledgment 

The  author  would  like  to  acknowledge  the  assistance  of 
Prof.  D.W.  Netzer  in  providing  information  which  was 
crucial  to  the  evaluation  of  the  analytical  results  obtained 
during  this  research. 

References 

•Warner,  C.F.  and  Netzer,  D.W.,  "Ac  Investigation  of  the  iTow 
Characteristics  of  Two-Phase  Flow  in  Converging-Diverging  Noz¬ 
zles,"  ASME  Paper  63-WA-192,  Winter  Annual  Meeting, 
Philadelphia,  Pa.,  Nov.  1963. 

1  Elliott,  D.Q.  and  Weinberg,  E.,  “Acceleration  of  Liquids  in  Two- 
Phase  Nozzles.”  NASA  Tech.  Kept.  32-967,  July  1968. 

’Clift,  R.,  Grace,  J.R.,  and  Weber,  M.E.,  Bubbles,  Drops  end 
Particles,  Academic  Press,  New  York,  1978. 

4 Wallis,  G.B.,  (Me-DimtHSkmtl  Two-Phase  Flow,  McGraw-Hill 
Book  Co.,  New  York,  1969. 


-59“ 


